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FORTRAN PROGRAM FOR CALCULATING VELOCITIES IN THE 
MERIDIONAL PLANE OF A TURBOMACHINE 
I - CENTRIFUGAL COMPRESSOR 
by Michael R. Vanco 
Lewis Research Center 

SUMMARY 

A FORTRAN IV computer program which calculates the velocities in the meridional 
plane of a centrifugal compressor is presented. This program will determine the ve- 
locities in the meridional plane of a backward -swept impeller, a radial impeller, and a 
vaned diffuser. The velocity gradient equation with the assumption of a hub-to-shroud 
mean stream surface is solved along arbitrary quasi -orthogonals in the meridional 
plane. These quasi -orthogonals are fixed straight lines. 

The input quantities for this program consist essentially of mass flow, rotational 
speed, number of blades, inlet total conditions, loss in relative total pressure, hub-to- 
shroud profile, mean blade shape, and a normal thickness table. The output yields 
meridional velocities, approximate blade surface velocities, streamline coordinates, 
blade shape coordinates, and stream -channel normal thickness in the meridional plane. 
Numerical examples are included to indicate the use of the program and the results ob- 
tained. 


INTRODUCTION 

Recently, increased interest has been shown in high-pressure -ratio backward -swept 
centrifugal impeller blades. Centrifugal compressors with backswept impeller blades 
have the potential of achieving higher efficiencies than those with radial impeller blades. 
Several methods are available for designing radial-bladed compressors, but limited 
work has been done on backward- swept impeller blades. Reference 1 gives the method 
and numerical techniques used to find the flow distribution in the meridional plane of a 
radial-flow turbine. This method solves the velocity gradient equations with the as- 
sumption of a hub-to-shroud mean stream surface. A set of arbitrary straight lines 



from hub to shroud is used instead of normals. These arbitrary straight lines are 
called quasi -orthogonals and they remain fixed regardless of any streamline change. 
This analysis, which has been used for radial-bladed centrifugal impellers, has now 
been programmed to include backward- swept centrifugal impeller blades. 

This report presents a computer program for calculating the velocities in the me- 
ridional plane of a centrifugal compressor. This program will determine the velocities 
in the meridional plane of a backward -swept impeller, a radial impeller, and a vaned 
diffuser, as well as approximate blade surface velocities. The output of this program 
is arranged in a form so that it can be used as input to programs used to calculate the 
blade -to-b lade loadings from references 2, 3, or 4. 

In this report, a description of the input and output and a FORTRAN IV computer 
program are presented. A brief description of the method of analysis and the computer 
program are given. Numerical examples are included to illustrate the use of the pro- 
gram and the results obtained. 


METHOD OF ANALYSIS 


Reference 1 presents the method and gives the numerical techniques used to find the 
flow distribution in the meridional plane of a radial-flow turbine. The general velocity 
gradient equation is derived along an arbitrary quasi -orthogonal in the meridional plane 
with the assumption of a hub-to-shroud mean stream surface. The equations derived in 
appendix B of reference 1 are 
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Figure 1. - Coordinate system and velocity components. 



The coordinate system and nomenclature are shown in figures 1 and 2. 

In this analysis, the total enthalpy at the inlet hj and the prerotation at the inlet X, 
r i^0i’ are assumec * constant. Therefore, equation (1) reduces to 

— = (a — + B — V + C — + D — (la) 

ds \ ds ds / ds ds 

Continuity must also be satisfied from hub to tip. The calculated mass flow across 
any fixed line from hub to tip must equal the specified mass flow. The mass flow is 
computed from 
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integrating from hub to tip along a quasi-orthogonal. 

The density is calculated from the isentropic flow equation with a correction for loss 
in total relative pressure. This equation is derived in reference 1. The density equation 
is 
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and the relative total pressure loss 


A P ,,= Pisen-P’ 


The change in the angular distance between blades A 0 is 
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where tangential thickness t^ is determined from 
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when the normal thickness t is specified. 

From figure 2, it can be seen that the velocity normal to the quasi -orthogonal is 

W n = W m cos(i// - a-) 


where from figure 1 


W m = w cos i 3 
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The flow angle 0 is determined from the mean stream surface, 9 = 0(m), for each 
streamline, between the blades. Therefore, 

sin a + ^ cos (12) 

where (d0/dm)j is the directional derivative along a streamline. 

The d6/d z and 30/3r in equation (9) refer to the mean blade shape. The (30/3z)j 
and (30/3r)^ in equations (la) and (12) refer to the mean stream surface between the 
blades. The mean stream surface is assumed to deviate from the mean blade shape at a 
radius r b for a centrifugal machine. An approximate equation for determining r b is 
given by reference 5, 



r b = r i e 


-0. 71(A0) 


(13) 


The equation for the mean stream surface when r > r b is 
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The boundary conditions used to obtain equation ( 14) were / 3 Q , the outlet flow angle; 0 b , 
the angular coordinate of the mean blade shape at r b ; and (d0/dm) b . Differentiating 
equation (14), we obtain 




— j(m - m b ) 2 




(15) 


It will be noted that equation (la) is in terms of (30/3z)j and (30/3r)^ and that, on 
the mean stream surface, 9 is a function of the meridional distance m, for each 
streamline. The relation between them is 



(16) 
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The preceding equations are solved with the specification of a mean blade shape. 

The mean blade shape can be specified by two methods. The first method of specifying 
the mean blade shape is specifying the angular coordinate of the mean blade shape 0 
constant along a quasi -orthogonal. Since the quasi-orthogonal is a fixed straight line, 
the mean blade shape is completely specified by specifying 0 as a function of the merid- 
ional distance m for the hub and shroud streamlines. Therefore, d0/dm is known, 
but the 89/dr and 30/3z have to be determined. If the directional derivative is taken 
in the m and s direction, then 
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and 
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With the specification of d0/ds = 0 and the geometry in figure 3, the following equations 
are obtained: 


30 _ cos xjy d 0 
3z cos(i \d - a) dm 

and 

30 _ sin ip d 0 
dr cos {\jy - a) dm 

This case is used for backswept centrifugal impeller blades. This case is also used for 
centrifugal diffusers, but equations (13) to (15) are not used because the mean blade 
shape is the same as the hub-to-shroud mean stream surface. 

The second method of specifying the mean blade shape is specifying 0 as a function 
of the axial distance z. This case is used for radial-element centrifugal impellers. 
Therefore, 30/3r = 0 and 
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Figure 3. - Relationship between coordin- 
ate systems. 


However, when slip occurs, that is, when the mean stream surface deviates from the 
mean blade shape, (d0/dm)j is known from equation (15). It is assumed that the mean 
stream surface deviates from the mean blade shape only in the radial direction. There- 
fore, dd/dz is known (mean blade shape), and 



( 21 ) 


The numerical techniques and procedures used for the solution of equations (la), 
(2), and (3) are given in reference 1. 


DESCRIPTION OF INPUT 

A description of the input for the FORTRAN IV computer program QUAC is given in 
this section. The input quantities consist essentially of mass flow, rotational speed, 
number of blades, specific -heat ratio, inlet total temperature and density, gas constant, 
loss in total relative pressure, hub-to-shroud profile, mean blade shape, and a normal 
thickness table. Since the program does not use any constants which depend on the sys- 
tem of units being used, any consistent set of units may be used. In the following input, 
each item has units specified in both the SI and U. S. customary systems. 


V 
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The input format is shown in table I. The first card is a title card and this card 
must be put in. The input variables are • 

MX number of quasi -orthogonals 

KMX number of streamlines 

MR number of r-values of TN in the thickness table 

MZ number of z-values of TN in the thickness table 

W rotational speed, rad/sec 

WT mass flow, kg/sec; slugs/sec 

XN number of full blades 

GAM specific-heat ratio 

AR gas constant, J/(kg)(K); (ft)(lbf)/(slug)(°R) 

TYPE integer; used as a code to indicate how arrays WA, Z, R, and DN are given 
initially; the integer values are 

0 These quantities will be calculated by the program. 

1 Quantities just computed for previous case will be used for next 

case. (Used only when more than one case is calculated on single 
computer run. ) 

MT number of z-coordinates in ZT array 

SRW integer that will cause the program to print out certain values; used for de- 

bugging purposes; the integer values are 

0 value when not debugging; usual case 

13 SPLINE 

16 SPLINT 

21 RUUT 

MXBL quasi -orthogonal number where blade starts 

TEMP inlet total temperature, Tl, K; °R 

2 2 

ALM inlet prerotation, \, m /sec; ft /sec 

RHO inlet total density, p!, kg/m^; slugs/ft^ 

2 2 

PLOSS loss in relative total pressure, Ap", N/m ; lb/ft 
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ANGR 


KSTH 


NPRT 


ITER 


streamline rotation angle, deg (The streamlines are rotated so that the slope 
of the program's cubic spline curve is not too large. Good results are ob- 
tained from the cubic spline if the absolute value of the slope is not greater 
than 1. Recommended angles are as follows: for an impeller, 45°; fora 
diffuser, 90°; and for an axial-flow compressor, 0°. ) 

determines the number of times the streamlines are smoothed for each itera- 
tion (For example, if KSTH = 0, no smoothing occurs. This is the usual 
case (KSTH = 0).) 

output control that determines which streamlines are printed out (For ex- 
ample, if NPRT = 1, every streamline is printed out; and if NPRT = 5, 
every fifth streamline is printed out. ) 

number of iterations to be performed after ERROR is less than TOLER or 
after ERROR has started to increase (If ITER = 0, data will be printed for 
every iteration; if ITER > 0, data will be printed only for the final itera- 
tion. Normally ITER = 1, but for a first-run set ITER = 0 and check the 
first few iterations to see if the data were put in properly. ) 


KD determines compressor type (For a backward -swept impeller, KD = 0; for 

a diffuser and an axial-flow compressor, KD = 1; for a radial element im- 
peller, KD = 2. ) 

SFACT blade multiplier to allow for splitter blades (For the case with no splitters, 
SFACT = 1. 0; and for the case with splitters, SFACT = 2.0.) 

ZSPLIT z-coordinate where splitter blade begins, m; ft (If there are no splitters, 

Z SPLIT > ZH(MX). ) 


BETO outlet flow angle, @ Q , deg 

CORFAC ratio of streamline correction used to calculated streamline correction 

(CORF AC affects the stability of the solution. If too large a value is used, 
the new streamlines are less smooth than the previous ones. If a computa- 
tion is based on this set of streamlines, the calculated streamline correc- 
tion becomes erratic. Therefore, it is important that the streamline cor- 
rection used give a smooth streamline for the next iteration. A value of 
0. 1 is recommended. ) 

SSN last quasi -orthogonal where smoothing is desired (For no smoothing, 

SSN = 0. ) 

ZS array of z-coordinates on shroud of hub-to-shroud profile located at quasi- 

orthogonal positions (see fig. 4), m; ft 
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RH 

THTA 

ZT 

TN 

XZ 

XR 
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« z 

Figure 4. - Input nomenclature. 


array of z -coordinates on hub of hub -to- shroud profile located at quasi- 
orthogonal positions (see fig. 4), m; ft 

array of r-coordinates on shroud corresponding to ZS (see fig. 4), m; ft 

array of r-coordinates on hub corresponding to ZH (see fig. 4), m; ft 

array of 6 -coordinates (mean blade shape), rad (When KD = 0 and KD = 1, 

6 is constant along a quasi -orthogonal and must correspond to the ZS, 

ZH, RS, and RH arrays. When KD = 2, 9 is a function of axial distance z 
and must correspond to the ZT array. ) 

array of z -coordinates corresponding to the THTA array, m; ft (Only used 
when KD = 2. ) 

array of thicknesses normal to the mean blade shape, t n , m; ft (This array 
has z-values of thickness going across and r-values of thickness going 
down the table. Values of thicknesses and corresponding z- and r- 
coordinates should extend beyond all boundaries of hub-to- shroud profile 
so that valid interpolation can be done in the program. ) 

array of z -coordinates for thickness table (TN), m; ft (The z -coordinates 
increase going across the table for a given r-coordinate. ) 

array of r-coordinates for thickness table (TN), m; ft (The r-coordinates 
increase going down the table for a given z -coordinate. ) 
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INSTRUCTIONS FOR PREPARING INPUT 


Theta Constant Along a Quasi-Orthogonal 

After the hub-to-shroud profile has been specified (fig. 5), the mean blade shape is 
determined. The angular coordinate of the mean blade shape 9 is specified as a func- 
tion of the meridional distance m for the hub and the shroud, as shown in figure 6. 
Values of 6 that are spaced to give good results from a cubic spline used in the program 
are selected. For a given value of 9 , the meridional distances are determined for the 
hub and shroud from figure 6. These meridional distances are then converted to the 


r 



► Z 

Figure 5. - Hub-to-shroud profile. 
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proper z- and r- coordinates. Therefore, the z- and r-coordinates for the end points 
of a quasi-orthogonal have been determined. These are the quantities 9 , r g , z g , r h , 
and z^ that are put in the program. The maximum number of quasi-orthogonal al- 
lowed is 21. 


Theta Not Constant Along a Quasi-Orthogonal 

This case is used for a radial impeller. The quasi-orthogonals are arbitrarily 
selected on the hub-to-shroud profile. They should be selected so that the program's 
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cubic spline curve will fit them smoothly. The mean blade shape is determined by spec- 
ifying 9 as a function of the axial distance z, as shown in the third numerical example 
(p. 18). MT is the number of 9 -values used. It should, also, be noted that KD = 2 for 
this case. 


Smoothing of Streamlines 

If the streamlines are not smooth, a smoothing routine can be used. KSTH is the 
number of times the streamlines are smoothed, and SSN is the last quasi-orthogonal 
where smoothing occurs. For an impeller, the streamline smoothing can take place 
only in the area shown in figure 7. It cannot take place in the other region because of 
the methods used. A recommended value for KSTH for smoothing is 4. 



I — 2 

Figure 7. - Streamline smoothing region for centrifugal impeller. 


Another method of smoothing the streamlines is to put quasi -or thogonals upstream 
of the impeller. The mean blade shape is extended into this region with the requirement 
of a negligible blade loading. These upstream quasi -orthogonals will allow a smoother 
transition into the impeller. For this case, MXBL is set equal to the quasi-orthogonal 
number where the blade starts. The first numerical example (p. 14) uses both these 
techniques. 


12 


DESCRIPTION OF OUTPUT 


An example of the output from the program is shown in table n. This output is in 
U. S. customary units. Each section of the output has been numbered to correspond to 
the following description: 

(1) The first output of the program is the input. 

(2) Output 2 gives the stagnation speed of sound at the inlet in meters per second 
(ft/sec); the radius at which the mean stream surface deviates from the mean blade 
shape (RB) in meters (ft); and a list of the number of iterations required to obtain a so- 
lution with the corresponding maximum streamline change in meters (ft). 

(3) Output 3 gives some of the important quantities used in the calculation procedure 
which are also useful for debugging purposes. This output is given for every streamline 
printed out. Streamline 1 is at the hub and streamline 21 is at the shroud. The number 
of streamlines printed out is controlled by the input parameter NPRT. Items listed are 


ALPHA 

angle between meridional streamline 

RC 

curvature of meridional streamline, 

SM 

meridional distance, m; ft 

BETA 

flow angle, j3, deg 

TT 

tangential blade thickness, m; ft 

SA 

A, eq. (2) 

SB 

C, eq. (2) 

SC 

B, eq. (2) 

SD 

D, eq. (2) 


(4) Output 4 gives the velocities and pressure for every streamline printed out. 
Items listed are 


Z 

R 

WA 

PRESS 

WTR 

WL 

TTREL 


z- coordinate, m; ft 
r-coordinate, m; ft 

relative velocity on mean stream surface, m/sec; ft/sec 
static pressure, N/m ; lb/ft 
suction-surface velocity, m/sec; ft/sec 
pressure -surface velocity, m/sec; ft/sec 
total relative temperature, K; °R 
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(5) Output 5 gives the stream -channel coordinates and the blade shape coordinates 
for the hub, mean, and shroud. Only the shroud information is shown here. This infor- 
mation is used to determine the blade-to-blade loading from reference 2, 3, or 4. The 
M ARRAY, R ARRAY, and the stream-channel normal thicknesses in the meridional 
plane are in meters (ft); and the THETA ARRAY, the angular coordinates of the blade 
shape, is in radians. 

STGR angular distance from center of trailing-edge circle of blade to center of 
leading-edge circle of blade, rad 

RI leading-edge radius, m; ft 

RO trailing-edge radius, m; ft 

For the case with splitters, the following additional output is given: 

MLER distance from leading edge of blade to leading edge of splitter, m; ft 

STGRS angular distance from center of trailing-edge circle of splitter to center of 
leading-edge circle of splitter, rad 

RI leading-edge radius of splitter, m; ft 

RO trailing-edge radius of splitter, m; ft 

BETAS flow angle at leading edge of splitter, deg 

(6) Output 6 gives the inlet flow angle for the hub, mean, and tip, in degrees. These 
angles are calculated inside the blade passage. 


NUMERICAL EXAMPLES 

To indicate the use of the program and the results obtained, three numerical ex- 
amples are given. The first example is a backward -swept centrifugal compressor rotor, 
the second is a centrifugal compressor diffuser, and the third is the input for a radial 
compressor. All examples are in U. S. customary units. 


Backward-Swept Centrifugal Compressor 

This compressor has a 6-to-l pressure ratio. The hub-to-shroud profile of the im- 
peller is shown in figure 8. The mean blade shape is given in figure 9, where 9 is 
specified as a function of the meridional distance m for the hub and shroud. The quasi- 
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Angular 




orthogonals shown in figure 8 depend on the mean blade shape in figure 9 because 9 is 
constant along a quasi -orthogonal . It will be noted that in this example three quasi- 
orthogonals were put upstream of the impeller. This was done to allow a smooth flow 
transition into the impeller because of the low inlet hub-to-tip radius ratio and the high 
rpm. Streamline smoothing was also used. MXBL was set equal to 4, SSN set equal to 
8. 0, and KSTH set equal to 4. The input for this case is given in table HI. The mean 
stream surface relative velocities are plotted in figure 10 for the hub, mean, and shroud 
streamlines. The velocity change near the impeller inlet was due to the blade blockage. 



The blade shape coordinates and the stream -channel normal thickness needed for calcu- 
lating the blade loading from reference 2, 3, or 4 are given in table IV for the mean 
streamline. These results are also obtained for the hub and shroud streamlines, but 
they are not shown here. 


Diffuser 

A flat-vaned diffuser for a centrifugal compressor was designed to have a linear 
static-pressure gradient from inlet to outlet. The meridional profile is shown in fig- 
ure 11. The angular coordinate of mean blade shape 9 is given as a function of the me- 
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Angular coordinate, 


I 



0 .01 .02 .03 

Axial distance, z, ft 


Figure 11. - Hub-to-shroud profile of compressor diffuser. 






ridional distance m in figure 12. The quasi -orthogonals shown in figure 11 depend on 
the mean blade shape in figure 12 because 9 is constant along a quasi-orthogonal . The 
input for this case is given in table V. The mean stream surface velocities and the ap- 
proximate blade surface velocities are plotted in figure 13 for the hub, mean, and shroud 
streamlines. The blade shape coordinates and the stream -channel normal thickness 
needed for calculating the blade loading from reference 2, 3, or 4 are given in table VI 
for the mean streamline. 

Radial Impeller 

This example is used to indicate the different input required. A hub-to-shroud pro- 
file is given in figure 14. The quasi -orthogonals for the profile shown are arbitrary and 
do not depend on the mean blade shape; that is, 9 is not constant along a quasi- 




Figure 13. - Blade loading diagram for hub and shroud of 
compressor diffuser. 
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orthogonal. The mean blade shape is put in as a function of the axial distance z, as 
shown in figure 15. Sample input is shown in table VII. The output obtained is the same 
as in the other examples. 


PROGRAM DESCRIPTION 
Main Program QUAC 

The main program QUAC contains all the equations given in the method of analysis 
and makes the majority of the calculations. It will be noted that K is used for the 
streamline number and I is used for the quasi -orthogonal number. QUAC calls the 
subroutines RUUT, SMOOTH, INTGRL, CONTIN, SPLDER, SPLINE, LININT, and 
SPLINT to perform various functions such as smoothing, finding roots, integration, in- 
terpolation, and use of a cubic spline curve to determine derivatives. These subroutines, 
excluding RUUT and SMOOTH, are described in reference 1. A brief description of each 
is given herein. 

The program variables for QUAC are 

A temporary storage 

AB temporary storage 
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AC 

temporary storage 

AD 

temporary storage 

AE 

meridional length from leading ec 

AL 

a 

alm 

see input 

AMLER 

MLER (see output) 

ANGR 

see input 

AR 

see input 

B 

temporary storage 

BA 

total weight flow between hub and 

BETA 

13 

BETAD 


BETAS 

see output 

BET AT 

% 

BETOH 

exit blade angle at hub 

BETOM 

exit blade angle at mean 

BE TOT 

exit blade angle at tip 

C 

temporary storage 

CAL 

COS O' 

CBETA 

COS |8 

Cl 

stagnation speed of sound at inlet 

CORF AC 

see input 

COSBD 

cos /3^ 

COSBT 

cos /3 t 

CP 

C P 

CURV 


DELBTA 

Pfh 

DELTA 

calculated streamline correction 


streamline 
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DENSTY 

DN 

DRDM 

DTD MB 

DTDMS 

DTDR 

DTDZ 

DWMDM 

DWTDM 

E 

ERROR 

ERROR1 

EXPON 

G 

GAM 

HR 

HZ 

I 

IND 

INF 

ITER 

K 

KD 

KMX 

KMXM1 

KSTH 

MR 

MT 

MX 


P 

distance along quasi-orthogonal from hub 

j 

(ro > + W sin 0)r A 6 

dm 

(dQ/ dm) b 

d0/dm at splitter leading edge 

de/dr 

de/dz 

dW /dm 
m' 

dW^/dm 

temporary storage 

maximum calculated streamline correction for present iteration 
ERROR from previous iteration 

l Ay - i) 

temporary storage 

y 

increment along quasi-orthogonal in r-direction 
increment along quasi-orthogonal in z-direction 
subscript to indicate number of quasi-orthogonal 
code number for use by subroutine CONTIN 
set equal to 1, when (ZH - ZS) = 0 
see output 

subscript used to indicate streamline number 

see input 

see input 

KMX - 1 

see input 

see input 

see input 

see input 
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MZ 

N 

Nl 

N2 

NPRT 

PLOSS 

PRS 

PSI 

R 

RB 

RC 

REXIT 

RH 

RHO 

RI 

RIS 

RO 

RS 

RSPLIT 

RUNO 

SA 

SAL 

SB 

SBETA 

SC 

SD 

SFACT 

SLA 

SM 


see input 
MXBL 
N + 1 
MX - 1 
see input 
see input 
P 
* 
r 

r b 

V r c 

average exit radius 
see input 

T 

Pi 

leading- edge radius 
leading- edge radius of splitter 
trailing-edge radius 
see input 

r-coordinate at leading edge of splitter 
run number 
A, eq. (2) 
sin a 

C, eq. (2) 
sin /3 

P, eq. (2) 

D, eq. (2) 
see input 

average distance between streamlines on a quasi-orthogonal 
distance from inlet along a meridional streamline 
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SMI meridional distance from first quasi-orthogonal to quasi- orthogonal that is be- 

fore point where stream surface deviates from blade surface 

SM2 meridional distance from first quasi-orthogonal to a quasi-orthogonal that is 

after point where stream surface deviates from blade surface 

SMF fractional meridional distance 

SMEXIT meridional distance from first quasi-orthogonal to trailing edge of blade 

SMRB meridional distance from first quasi-orthogonal to point where mean stream 

surface deviates from mean blade shape 

SRW see input 

SSN see input 

STGR see output 

T t (interpolated value) 

TANBB tan ^ 

TANS tan /3 , at leading edge of splitter 

D 

TEMP T[ 

THTAB 0 b 

THTAF 0 f 

THTAS 0 

D 

THH mean blade shape 0 -coordinate at hub 

THHC temporary storage 

THHl blade shape, 0 -coordinate at hub on surface 1 
THH2 0 -coordinate at hub on blade surface 2 
THM 0 -coordinate of mean blade shape at mean 

THMC temporary storage 

THMl 0- coordinate at mean on blade surface 1 

THM2 0- coordinate at mean on blade surface 2 

THS 0- coordinate of mean blade shape at shroud 

THSC temporary storage 

THSI 0 -coordinate at shroud on blade surface 1 

THS2 0-coordinate at shroud on blade surface 2 
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TN see input 

TOLER iteration tolerance 

TSPLIT normal blade thickness at leading edge of splitter 
TPP1P T"/T! 

TTREL see output 

TT t e 

TYPE see input 

TIP T/T! 

W u 

WA W 

WAS W*, eq. (13) of ref. 1 

WASS W**, eq. (13) of ref. 1 
WT total mass flow 

WTFL calculated total mass between hub and streamline 

WTHRU W„ 

n 

WTR W t , eq. (10) of ref. 1 (suction- surface velocity) 

WL pres sure -surface velocity 

XN see input 

XR see input 

XZ see input 

YA average weight flow per unit length crossing a quasi -orthogonal 

YH temporary storage 

YM temporary storage 

YS temporary storage 

Z z 

ZEXIT average z-coordinate at exit 

ZH see input 

ZS see input 

Z SPLIT see input 

ZT see input 
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Subroutine RUUT 


Subroutine RUUT finds the root between two given points. It is used to find the me- 
ridional distance where the mean stream surface deviates from the mean blade shape 
when the radius at which this occurs is given. If the root cannot be found within the tol- 
erance, a message is printed out and the input arguments are listed. If-there is trouble 
in finding a root, set SRW = 21 in the input and all the input to the subroutine will be 
printed out. 

The calling sequence for RUUT is 

CALL RUUT(SM1, SM2, RB, SMRB, SM(1, K), R(l, K), MX) 


where 


SMI 

SM2 

RB 

SMRB 
SM(1, K) 
RU,K) 
MX 


meridional distance of quasi-orthogonal before point desired (input) 

meridional distance of quasi-orthogonal after point desired (input) 

radius at point desired (input) 

desired meridional distance (output) 

array of m- coordinates (input) 

array of r-coordinates (input) 

number of r-coordinates (input) 


Subroutine SMOOTH 

Subroutine SMOOTH smoothes the streamlines to obtain abetter numerical solution. 
It uses the hub streamline as the base streamline for the smoothing operation. 

The slopes of the quasi- or thogonals are 


m 


I 



and the streamline slopes are 


m 


K " 


?i + i - yi-i 

x i+i ■ x i-i 


( 22 ) 


(23) 
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where K is the streamline number and I is the quasi -orthogonal number. The x-y 
coordinates of an intersection can now be determined. From analytical geometry, 


(yj_l - m K x I-l) ' ( y h T “ m K x h T ) 

x, = 

' L 'l mj - m^. 

The smoothed x-coordinate is 



where D is the smoothing factor. The smoothed y-coordinate is 

y, = 

When mj = 0, the following equations are used: 

(yi) l = m K< x i- x i-i> + yi-i 



and 

( X l) x = X I 

When m^ = 0, 

( y l) I = m l( X I-l" X h I ) + y h I 



( 24 ) 


(25) 


(26) 


(27) 


(28) 


(29) 


( 30 ) 
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and 



( 31 ) 


The value of D is 2 for all quasi- orthogonals except for the last three, where 
smoothing occurs. The values of D for these three are 2. 6667, 4.0, and 8.0, respec- 
tively. This was done so that there would not be any discontinuities when only certain 
sections of the streamlines are smoothed. 

The calling sequence for SMOOTH is 

CALL SMOOTH(Z(l, K), R(1,K), ZH, RH, AB, SSN, INF) 

where 

Z(l, K) z-coordinate of streamline 

R(l, K) r-coordinate of streamline 

ZH z-coordinate of hub streamline 

RH r-coordinate of hub streamline 

AB slope of quasi-orthogonals, mj 

SSN last quasi-orthogonal where smoothing is desired 

INF indicator for quasi-orthogonals with a slope of infinity 

The program variables are 


D 

SLOPE 
SLOPE 1 

X 
XH 

XI 
Y 
YH 
Y1 


smoothing factor 
slope of quasi-orthogonals, nij 
streamline slope, m^ 
z-coordinate of streamlines 
z-coordinate of hub streamline 
z-coordinate of smoothed streamline 
r-coordinate of streamlines 
r-coordinate of hub streamline 
r-coordinate of smoothed streamline 
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Other Subroutines 


Subroutines INTGRL, CONTEST, SPLDER, SPLINE, LININT, and SPLINT are de- 
scribed in reference 1. INTGRL is used for numerical integration. CONTEST is used to 
determine the hub velocity for the next continuity iteration. SPLDER is used to deter- 
mine the values of the derivatives at the specified interpolated points. SPLINE is used 
to determine the first and second derivatives. If there is a problem with the SPLINE 
subroutine, set SRW = 13 in the input and the input and output of the SPLINE subroutine 
will be printed out. LININT is used to determine the interpolated values of the normal 
blade thickness from the given thickness table. SPLINT is used for interpolation. The 
input and output data for SPLINT will be printed out if SRW = 16. 


PROGRAM LISTING 


$ I RF TC GUaC DECK 

C CaLCULATI ON OF VELOCITY AND PRESSURE D I STR I BUT I ON IN A CFNTR I FUSAL COMPRESSOR 
C PY USE OF GUAS I -QRT HOGON AL S 

C 
C 

COMMON SRW 

01 MENS I ON AL(2l,21),BFTA(21,21),CAL(21,21),C8ETA(21,21),INF(21», 
1CURV(21,22) ,DN( 2 1,21), PRS< 21,21 >,R( 21,2 l >,Z( 21,21 >,SM( 21,21 > , 

2SAI21 ,21 ) ,SR(21, 21 ) .SCI 21, 21 I, SOI 21,21 ) ,SAL (21,21) »SRETA(21 ,21 ) , 

3TNI21.21 ) ,TT(21, 21 ) ,WA( 21,21 ),WTR{ 21,21) , TTREL ( 2 1 , 2 1 ) , WL (21 ,21) 

DIMENSION AR( 22) , AC( 2?) , A0( 22), BA( 2 1 ) , DELS T A ( 2 1 ) ,DRDM(21 ) , AM??) , 

1 YM (21), 0TDMI21 ) ,DWMDM( 21 ) ,0WTDM( 21 ),RH( 21 ) ,RS( 21 >,ZH(21),ZS(2l), 
2THTA(21),WTFL(21),XR<21).XT<21),XZ(21),RETAI(3),AA(3>,THTAF(21> 

DIMENSION THH(2l ) ,T HM ( 2 1 ) , THS ( 2 1 ) , THH 1 ( 2 1 ) , THH2 (21 ) , THM1 (21 ) , 

1THM2I21 ) ,THS1 <21 ),THS2( 2) ) ,DTDZ (21 ),DTDR(21 ) ,ZT{ 21) , 

2YA < 21 ) , YM 21 ) ,YS (21),Tim,TO(3t 
INTEGER R UNO, TYPE, S RW , HUB , SHROUD 

RUNC=0 

10 READ (5,1001) 

WRITE<6,1049) 

WRITE (6 , 1001 ) 

READ t 5, 1010) MX , KMX, MR ,MZ , W , WT , XN , GAM , AR 

ITNC = 1 

RUN0= RUNC + 1 

WRITE (6,1020) RUNO 

WRITE(6,1007) 

WRITE (6,1011 ) MX, KMX, MR, MZ , W , WT , XN , C AMI , AR 

READ (5,1010) TYPE *. MT , SRW , M XRL , TEMP , AL M , RHO , P L 0 SS , ANGR 

WR I TE ( 6 , 1008 ) 

WRITE (6, 1011) TYPE, MT, SRW, Ml XRL, TEMP, AL M , RHO , PLOSS , ANGR 

RE AC (5,1010) KSTH,NPRT , ITER»KD , SF ACT, Z SPL I T, BETO ,CORFAC,SSN 

WRITE (6, 1009) 

WR I TE ( S , 10 1 1 ) KSTH.NPRT, ITER.KD , sfact, zspl I T,BET0 ,CORFAC,SSM 
ITER1 = ITER 

RE A C ( 5 , 1C 30 ) ( ZS ( I ) , 1 = 1, MX) 

WR I TE ( 6 , 1 029 ) 

WR I TE ( 6 , 1028 ) (ZS ( I ) , 1=1 ,MX ) 

REAC(5, 1030 ) (ZH( I), 1=1, MX) 
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WRITF(6,1031) 

WRITE (6, 1028) ( ZH ( I ),1=1,MX) 

RE AC (5, 1030 MRS! I) , 1=1, MX ) 

WRITEI6, 1032) 

WRI TE( 6 ,1028) IRS (I) , 1 = 1 , MX ) 

READ (5, 1C 30) ( RH ( I) , 1=1, MX) 

WRITE(6»1033) 

WRITE (6, 1028) t RH III , 1 = 1, MX ) 

IF(TYPF.NE.O) GO TO 145 
IFIRSIl ) .EC.RH(l) > GO TO 20 

WA ( 1 , 1 ) = WT/RHO/l (RS< 1)**2-RH( 1)**2)*3.14) 

GO TC 21 

20 WA(1 ,1) =WT/RHO/( ZH( 1)-ZS( 1) 1/3.14/ (RSI 1)+RH(1) ) 

21 DO 30 1 = 1, MX 

DN ( I , KMX ) = SQR T ( ( ZS( I )-ZH( I ) )**2+(RS( I )-RH( I ) )**2) 
DO 30 K=1 , KMX 

DN ( I ,K)=FLCAT(K-1)/FL0AT(KMX-1)*DN( I, KMX) 

WA( I ,K)=WA(l,l) 

Z ( I ,K)=DN( I ,K)/DN( I,KMX)*( ZS ( I )-ZH( I ) !+ZH( I ) 

3C R ( I ,K) = DN ( I ,K )/DN( I, KMX )*( RSI D-RHI I I )+RH( I ) 

IF ( KD.EQ.2 ) GO TO 50 
READ I 5,1030)!THTA( I ), 1=1, MX ) 

WR I TE I ft , 1034 ) 

WRITE (6,1028)(THTA( I), 1 = 1, MX ) 

WRITE(6,1035) 

GO TC 51 

5C REAO (5,1030)(THTAI I), 1=1, MT ) 

WR I TE l 6 , 1034 ) 

WRITE (6,1028 ) ITHTAI I ), 1=1, MT ) 

RE AD I 5, 10 30 ) ( ZT I I ) , I = 1, MT ) 

WR I T E I 6 , 1039 ) 

WRITE 1 6, 1028) I ZT (I) , I=sl , MT ) 

51 WR I TE I 6 , 1 036 ) 

DO 60 K = 1 , MR 

RE AC (5,1030)(TN( I ,K) , 1=1, MZ) 

6C WRITE (6,1028) (TNI I»K), 1 = 1, MZ) 

RE AC ( 5,1030) (XZ( I),I=1,MZ) 

WRITE(6,1037) 

WRITE ( 6, 1028 ) (X7( I ) , 1 = 1, MZ ) 

READ (5 , 1030 ) (XRI I ) , 1=1, MR ) 

WRITE (ft , 1038 ) 

WRITE ( 6, 1028 ) ( X R { I ) , 1 = 1, MR) 

END CF INPUT STATEMENTS 

INITIALIZE , CALCULATE CONSTANTS 

WTCLER = WT / 1000 00 • 

TCLER = (RS(1 )-RH(l) 1/5000. 

IFIRSIl). EC. RH<1>) TOLER= ( Z H l 1 ) — Z S ( 1) 1/5000. 

DC 110 K= 1 , KMX 
1IC SM(1,K)=C. 

BA(l)=iO. 

DC 120 K= 2 , KMX 

1 2 C B A ( K ) = FLCAT(K-1)*WT/FL0AT(KMX -1) 

DC 1 30 1=1, MX 
13C DN ( I , 1 ) =C • 

ANGR = ANGR/57. 29577 
145 CONTINUE 

Cl = SORT < C-AM#AR*TEMP ) 

WR I T E ( 6 , 1049 ) 

WRITE (6,1050) Cl 


KMXM1 = KMX-1 
CP= AR*GAK/ < GAM-1 . ) 

FXPCN = 1. /(GAM-1.) 

BETC = BETO /57. 29577 
ZEXIT = <ZS (MX)+ZH(MX ) 1/2. 

REXIT = (RS(MX)+RH(MX) )/2. 

IF ( KD.EQ. 1 ) GO TO 149 

CALL LININT (ZEXIT , REXIT , X Z , XR , TN , 2 1 , 2 1 , T ) 

RB = REXIT *EXP(-.71*(2.*3.14159/(XN*SFACT)-T/REXIT )) 

WRITE (6,1027) RB 
149 ERRCR = 1 OC 0C0 • 

C 

C BEGINNING OF LOOP FOR ITERATIONS 

15C IF( ITER.EQ.O) WRITE (6,1060) ITNO 
C 

I F ( ITER.EQ.O) WRITE (6, 1070) 

ERRCR 1= EP RCR 
ERRCR=0 . 

C 

C START CALCULATION OF PARAMETERS 
C 

DC 180 K= 1 * KMX 
DC 180 1 = 2, MX 

SM ( I , K ) = S M ( I — 1 , K ) + SQR T ( ( Z ( I , K )-Z ( I — 1 , K > )**2+(R( I , K ) -R ( I — 1 , K ) > ** 
1 2) 

180 CONTINUE 

DC 230 K = 1 , K M X 
DC 160 1=1, MX 

A B ( I ) = Z ( I ,K )*C0S< ANGR ) + R( I * K ) *S IN( A NCR ) 

160 AC ( I ) = R ( I,K)*COS(ANGR) - Z ( I , K ) *S I N ( ANGR ) 

CALL SPLINE ( AR , AC , MX , AL ( 1 , K ) , CURV ( 1 , K ) ) 

DC 170 1=1, MX 

CURV ( I,K )=CURV< I »K)/(1.+AL(!»K)**2)**1.5 
A L ( I , K ) = ATAN(AL( I , K > )+ANGR 
CAL ( I »K ) = COS( AL( I , K ) ) 

17C SAL ( I ,K ) = S IN( AL ( I , K ) ) 

IF ( KD.EQ. 2) GC TO 171 

CALL SPLINE ( SMI 1 ,K ) ,THTA , MX , DTDM , AC ) 

GC TC 172 

171 CALL SPLCER(ZT,THTA,MT,Z(1,K),MX,DTDZ) 

172 DC 2C4 I = 1 , MX 
T = 0. 

THTAF ( I ) = T H T A ( I ) 

I F < I.GE.MXRL) CALL L IN I NT( Z ! I , K ) , R ( I , K ) , XZ , XR , TN , 21 , 2 1 , T ) 

IF (ZS(I).GE.ZH(I)) C-OTO 202 

PSI = ATAM <RS( I )-RH( I ) )/( ZS( I )-ZH( I ) J ) + l. 5708 
GC TC 202 

202 PSI = ATA N< ( ZH( I )-ZS( I ) )/( RS ( I )-RH( I ) ) > 

203 IF ( KD.EQ. 2 ) DTDM(I) = DTCZ ( I ) *C AL ( I , K ) 

IF ( KD.EQ. 2 ) CTDR ( I ) = 0.0 

IF (KC.NE. 2 ) DTDZ(I) = COS ( PS I )/ COS ( PS I - AL ( I , K ) ) *DTDM (I) 

IF (KC.NE. 2 ) DTCR(I) = S IN ( PS I ) /COS ( PS I - AL ( I , K ) ) *DTDM ( I ) 

2C4 T T ( I , K ) =T*SQRT( 1.0 + R( I , K ) ** 2* ( DTDR ( I )**2+0T0Z( I )**2) ) 

IF (KC.EC.l) GO TO 207 
DC 205 I = 1 , M X 

IF (R(I,K).GT. RR) GO TO 206 

205 CONTINUE 

206 SMI = SM ( I - 1 , K ) 

SM2 = S M ( I , K > 

CALL RUUT ( SM1,SM2,RP,SMRB,SM( 1 ,K ) ,R( 1,K ) ,MX) 

IF (KC .EG. 2 ) CALL SPLINT (ZT,THTA»MT*Z(1,K), MX, THTAF) 

CALL SPLINT (SM( 1»K), THTAF, MX, SMRB,1,THTAB) 


noon 


CALL SPLCER (SM( 1 ,K ) , THT AF , MX , SMRB , 1 , DTDMR ) 

TANFB = RB*OTDMB 
SMFXIT = SM(MX,K) 

DC 2C1 I = 1,MX 

IF <R ( I ,K ) .LT.RB ) GO TO 201 

THTAF(I) = THTAB + ( S M< I , K ) - SMR B ) **3* ( T AN ( BETO ) ZRFX I T- TANBB /RB ) / 

1 I3..C* CSMEX IT-SMRB ) **2 ) + ( SH { I , K )-SMRB > * TANRB/RB 
DTCM(I) = + I SM( I ,K )-SMRB >**2*( TAN(RETO) /REX I T-TANBB/RB ) / 

1 ( ( SME X I T- SMRB )**2 I+TANBB/RB 
IF ( SAL(I.K) .EQ. 0.0 ) GO TO 200 
DTCR(I) = ( DTDM< I )-DTDZ< I)*CAL( I ,K )) /SAL ( I ,K) 

GC TC 201 

200 DTCR(I) =0.0 

201 CONTINUE 

207 DC 220 1=1, MX 

BETAII.K) = ATAMIRI I,K)*DTDMIII ) 

SBETA(I,K) = SIN ( BET A( I ,K ) ) 

CBETAII.K) = COS (BETA! I , K ) ) 

Afi ( I ) = VI A ( I ,K)*CBETA( I ,.K ) 

22C AC ( I )=WA(I»K)*SBETA( I,K) 

CALL SPLINE(SM(1,K),AB, MX, DWMDM, AD ) 

CALL SPLINE(SM(1,K) , AC , MX , OWTQM , AD ) 

IF( ( ITER. LE.O). AND. < MOO < K- 1 , NPRT >. EQ . 0 ) ) WRITE (6,1C80) K 
DC 23C 1=1, MX 

S A ( I , K ) = CRETA( l ,K ) **2*CAL < I,K )*CURV( I ,K l-SBETA ( I ,K)**2/R( I ,K) 

1 + SAL (I , K1*CBETA< I »K) *SRETA( l ,K) *DTDR ( I ) 

SB ( I , K ) = SAL(I,K)*CBETA(I,K) * DWMDM ( I ) - 2 . 0*W*SB ETA ( I , K ) +DTDRIII 
1*R(I,K)*CBETA(I,K )*< CWTOMI I I +2 . * W* S AL ( I , K ) ) 

SC ( I , K ) = -CBETA ( I ,K)**2*SAL ( I ,K )*CURV( I,K ) 

1+SAL(I»K)*CRETAI I,K)*SRETA( I,K )*DTD7.( I > 

SD(I.K) = CAL ( I, K )*CBETA( I,K ) *DWMDM II) +DTDZII) 

1*R( I ,K) *CRETA< I , K ) * ( DWT OM ( I ) +2.*W*SAL< l * K ) ) 

IF ( ( ITER.GT.O J.OR. ( MODI K-l , NPRT ) .NE.O) ) GO TO 230 
A= AU I,K)*57.29577 
8= S M ( I , K ) 

E= TT(I.K) 

G = R E T A ( I ,K)*57. 29377 

WRITF (6,1090) A,CURV( I,K),B,G,E» S A ( I , K ) , SR I I , K ) , SC ( I , K ) , SD ( I , K ) 
2 3 C CONTINUE 

END OF LOOP - PARAMETER CALCULATION 

CALCULATE BLADE SURFACE VELOCITIES (AFTER CONVERGENCE) 

IF( ITER. NE.O) GO TO 260 
DC 25C K= 1 , KMX 

CALL SPLINE (SM< 1»K) ,TT< 1,K) , MX , DEL BT A , AC ) 

A = X N 

DC 2 AO 1=1, MX 

A B 1 I ) = ( R( I ,K)*W + WA( I , K ) *SBET A ( I ,K ) )* ( 6.2831 86*R ( I ,K) / A-TTI I ,K) ) 
CALL SPLINE ( SM < l ,K ) , AR , MX , ORDM , AC ) 

IF (SFACT.LE. 1.0) GC TO 2A5 
A = SFACT *XN 
DC 2AA 1=1, MX 

AB ( I ) =( R ( I ,K) *W+WA( I,K)*SBETA(I,K) ) * I 6 . 28 3 1 86*R < I , K ) / A-TTII ,K) ) 
CALL SPLINE ( SM ( 1 , K ) , AR , MX , AD , AC ) 

DC 250 1=1, MX 

RETAC = EETA( I,K)-DELBTA( I )/2. 

BE T A T ' = BET AO + DE L BT A ( I ) 

CCSRC = CCS(BETAC) 

CCSPT = COS(BETAT) 

I F ( Z ( I ,K ) .GT. ZSPL IT ) DRDM(I) = ADI I ) 


2AC 

2AA 
2 A 5 
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WTR ( I ,K)=CCSRD*COSBT/(COSBD+COSBT)*{ 2.*WA( I ,K ) /C OSBO + R ( I , K ) * W* 

1 (BET AC- BET AT) /CB ETA I UK ) **2+CRDM ( I ) ) 

WL ( I t K ) = 2.0*WA ( I,K )-WTR< I,K) 

25C CONTINUE 

ENC OF BLADE SURFACE VELOCITY CALCULATIONS 
START CALCULATION OF WEIGHT FLOW VS. DISTANCE FROM HUB 

26 C DC 37C 1=1, NX 
I N C = 1 

DO 270 K= 1 , KMX 
27C AC(KI=DMItK) 

GC TC 29C 

28C WA(U1) = .5*WA(I, 1) 

29C DO 3C0 K = 2 » K N X 
J- K- 1 

HR=R ( I,K)-ft( I,J) 

HZ=Z < UKI-Z(UJ) 

WAS = WA( I , J)*( l.O+SAt I, J)*HR+SC( I* J )*HZ) + SB ( I , J ) *HR+SD ( I , J ) #HZ 
WASS =i WA (I , J)+W AS + ISAI I,K )*HR+SC( I,K) *HZ ) +SB( I,K )*HP+SD( I ,K)*HZ 
3CC WA(I,K)=(WAS+WASS)/2. 

3 1 C DO 3 AC K= 1, KNX 

T 1 P= 1.-<WA ( UK) **2+2.*W*ALM-(W*R( I,K) ) **2 > /2 . /C P/TEMP 
IFfTIP.LT. .0) GO TO 280 

TPP1P=; 1.- < 2.*W*ALM-(W*R( I,K) ) ** 2 ) / 2 . /C P /TEMP 

TTREL ( I ,K ) = TPP 1P + TEMP 
SNF = O.C 

I F ( UGE.NXBL) SMF= ( SM ( I, K )-SM ( MXBL , K ) ) / ( SM ( MX , K ) -SM ( MXBL , K ) ) 
DENSTy=T1P**EXP0N*RH0-( T1P/TPP1P)**EXP0N*PL0SS/AR/TPP1P/TEMP*SMF 
PR S < UK UOENSTY* AR*T1P*TEMP 
I F ( ZS (I) .GE.ZH< I )) GC TO 320 

PSI = ATANI (RSI I )-RH( I) )/(ZS< I )-ZH( I ) ) >+1.5708 
GC TC 33C 

32C PSI=ATAN((ZH(I)-ZS<I))./(RSm-RH(I)>) 

33 C WTHRU=WA( I ,K)*CBETA( I, K ) *COS ( PS I -AL ( I,K) ) 

A = XN 

IF (Z I UKJ.GT.ZSPUT) a=sfact*xn 
C = 6.282 186*RU ,K)-A*TT( I,K > 

34 C AD(K)=;DENSTY*WTHRU*C 

CALL INTCRL(AC(1UA0(1) ,KMX,WTFLI 1 ) ) 

Y A ( I ) = VTFUKMX )/DN( I, KMX ) 

YH< I ) = £0(1) 

KN = (K N > + 1 )/2 
YNtI) = AD(KN) 

YS ( I) = AD(KNX) 

IF (flPSIVT-WTFL(KMX) J.LE.WTOLER) GO TO 350 
CALL CONTIN ( WA ( I, 1 ) , WT FL I KMX ) , IND , I, WT ) 

IF (IND.NE.6) GO TO 290 
35C CALL SPLINT < WTF L , AC , K'M X , 8 A , KMX , AB ) 

DC 3ft 0 K= 1 , K NX 
DELTA=ABS ( A8(K)-CN( I » K! ) ) 

DN ( I , K ) = ( 1 .-CORF AC ) *CN ( I ,K ) +COR F AC*AB ( K ) 

36 C IF ( DELTA. GT. ERROR ) ERROR = OELT A 
37C CONTINUE 
C 

C END CF LOOP - WEIGHT FLOW CALCULATION 
C CALCULATE STREAMLINE COORDINATES FOR NEXT ITERATION 
C 

DC 380 K=2,KMXM1 
DC 3 SC 1 = 1, NX 

Z ( I,K ) = CN ( I,K)/ DNI I,KMX)*(ZS( I )-ZH< I)) «-ZH( I ) 
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38 C R(I,K) = DMI,K)/DN(I,KMX)*(RS(l )-RH( I ) )+RH( I ) 

IF (KSTH.EC.O) GC TD 383 
DC 381 1=1 , W X 
I N F ( I ) = 0 

IF(ZS(II.EC.ZH(I)I GC TO 3805 
APII) = (RS(II-Rnmi/(ZSUI-ZH(II) 

GC TC 381 
3805 INF(I) = 1 

381 CCNTIMJE 

DC 382 K=2 * K M X M, 1 
DC 382 J= 1 * KSTH 

382 CALL SMCCTH I Z ( 1 , K ) , R ( 1 . K 1 , Z F , RH , AB , SSN , I NF 1 

383 IF ( ( FRROP .GE.ERR0R1 ) .OR. I ERROR .L E. TOL ER > ) ITER = I TER-1 
I F ( ITER.CT.O) GO TO 410 

WRITE (6,1100) 

DC A C 0 K=1,KMX,NPRT 
WRITE (6,1080) K 
DC 390 1=1, MX 

AB I I ) = Z ( I,K)*CCS( ANGR) + R < I , K ) *S IN ( ANGR ) 

390 AC(I) = R ( I ,K)*COS( ANGR) - Z ( I , K ) *S I N ( ANGR ) 

CALL SPLINE ( AB, AC , MX ,.AD , CURV ( 1 , K ) ) 

DC 4C0 1=1, MX 

CURVI I , K ) = CURV ( I » W ) / ( 1 • + AD ( I ) **21 **1.5 

B= 7 ( I , K ) 

D= R ( I , K ) 

4CC write (6,1110) R,D,WA(I,K),PRS(I,K),WTR(I,K),WL(I,K),TTREL(I,K) 
WRITE (6,1130) 

41 C A= ERRCR 

WRITE (6,1120) I TNO , A 
I TN C = I T NC + 1 

IF ( ITER.GE.O) GO TO 150 
N = MXBL 
CC A 1 9 J = 1 , 3 
K = 1 

IF (J.EC.2) K = (KMX+ll/2 
IF (J.EQ.3) K = KMX 
IF (KC .EC. 2 ) GC TO 417 

CALL SPLINE (SM( 1,K),THTA ,MX,DTDM,AC) 

GC TC 41E 

417 CALL S°LCER(ZT,TFTA,MT,Z(1,K),MX,DTDZ) 

D T C M ( MX ) = CAL(MX,K)*CTOZ(MX) 

418 IF (J.FC.l ) BETCH = AT AN ( R ( MX , K ) *DTDM ( MX ) ) 

IF (J.EC.2 ) BET CM = AT AN ( R ( MX , K ) *DTDM ( MX ) ) 

IF (J.EG.3 ) BETCT = AT AN ( R ( MX , K ) *DTDM ( MX ) ) 

CALL LININT ( Z ( M X , K ) , R ( MX , K ) , X Z , XR , TN , 2 1 , 2 1 , TO ( J ) ) 

419 CALL LININT <Z( N,K),R( N , K ) , X Z , XR , TN , 2 1 , 2 1 , T I ( J ) ) 

K = (KMX+l)/2 

DC 440 1 = 1 , MX 

SLA = DM I ,KMX)/FL0AT(KMX-1) 

IF ( 7S( I ) .GE.ZHl I ) V GC TO 420 

PSI = ATAN( (RS( I )-RH( I) )/( ZS ( I )-ZH( I ) ) J + 1.5708 
GC TC 43C 

42C PSI = ATAN( (ZH( I ) — 7S( I) )/(RS( I )-RH( I ) ) ) 

430 A B ( I ) = V A( I ) *SL A*COS(PSI-AL ( I, 1 ) )/YH(I) 

AC(I) = V A ( I )*SLA*COS(PSI-AL( I , K ) )/YM( I ) 

440 AC < I ) = YA( I)*SLA*COS(PSI-AL(I»KMX) )/YS( I) 

IF ( KD .EC. 2 ) GO TO 442 

DC 441 1=1 , MX 

THF ( I ) = THT A ( I ) 

THM { I ) = THTA ( I ) 

441 THS ( I ) = TETA (I) 

GC TC 443 


442 


CALL SPLINT ( ZT , THT A f MT , 1 { 1 , 1 },MX,THH) 

CALL SPLINT < ZT , THT A , MT , Z ( 1 , K ),MX,THM) 

CALL SPLINT ( Z T , THT A , MT , Z I 1 , KMX ) , MX , THS ) 

443 R I = T I (1) /2. 

THHC= THHINI+RI *TAN(BETA(N, 1) ) /R( N * 1 ) 

rc = ram / 2 . 

THH1(MX)=THH(MX)-R0*TAN(BET0H)/R(MX,1 ) -THHC 
THH2(MX) = THH<M.X)-R0*TAN(BET0H)/R(MX,1 ) -THHC 
R I = TI (2 ) /2. 

THMC = TEM (N) +RI *TAN( BET A I N,K ) ) /R( N, K ) 

RC = TO ( 2 )/2. 

THM (MX )=THM( MX) -RO*TAN( BETOM) /R(MX,K ) -THMC 
TH M 2 ( MX ) = TH M ( MX ) -RO*T ANI RETOM)/R(MX,K ) -THMC 
R I = TI (2 I / 2 - 

THSC= THS(NH-Rl*TAN( RETAIN, KMX ))/R(N,KMX) 

RC = TOI2I/2. 

THS1 ( MX ) = TH S ( M X) -RO*T ANI BETOTl/R ( MX, KMX) -THSC 

THS2(MX)=THS(MX)-RO*TAN(RETOT)/R(MX,KMX) -THSC 

CC 449 I = 1 , MX RL 

THHI(I) = 0.0 

THM 1 I I ) = 0.0 

THS I ( I ) = 0.0 

THH2II) = 0.0 

THM2II) = 0.0 

449 THS2II) = 0.0 
N I = N+l 

N2 = MX - 1 
DC 450 I = N 1 , N 2 

THHKII = THHIIH-: TT(I *l)/2./RC 1,1 )-THHC 

TH M 1 ( I ) = T H M ( I ) +' TT(I,K)/2./R(I,K> -THMC 

THS 1 ( I ) = THSIII+ TT(I,KMX )/2./R(I,KMX ) -THSC 

THH2II) = THHIII- TTI I , 1 ) /2 . /R 1 I , 1 ) -THHC 

THM2II) = THMIII- TT < I ,K ) / 2./R ( I , K ) -THMC 

450 THS2II) = THS I I )— TTII,KMX >/2./R(I,KMX > -THSC 

WRI TE (6, 1200) 

WRITEI6.1239) 

WRITE (6,1251 ) 

DC 451 1=1, MX 

451 AE(I) = SMII.1 )-SM(MXBL,l ) 

WRITE (6, 1230 MAE (II , T=.l,MX) 

WP ITE (6 , 1249) 

WR I TE ( 6 , 1 230 ) ( R(I,1 ),I=1,MX) 

WR I TE (6 , 1 240 ) 

WRITEI6.12 30) ( A B ( I ) , I = 1 , MX ) 

WRITE(6,1250) 

WR l TE ( 6 , 1 25 1 ) 

WR ITE (6 , 1230) I AE (I), 1 = 1, MX ) 

WRITEI6.1252) 

WR I TE (6 , 1230) ITHH1 (I ) ,-1 = 1,MX ) 

WRITE (6,1253) 

WRITE (6, 1230) (THH2 ( I),. 1 = 1, MX) 

R I = TI ( 1 > /2. 

RC = TQ(1 ) /2. 

ST GR = TEHl(MX) 

WRITE (6, 1254) ST GR , R I ,RO 

IF ( Z ( M X , 1 ) .LT. ZSPLIT) GO TO 453 

CALL SPLINT (Z(l,l l*SM(l,l ) , MX , Z SPL I T, 1 , AML ER ) 

CALL SPLINT ( SM < 1,1 ),R(1,1 ) , MX , AML ER , 1 , RSPL I T ) 

CAUL SPLINT (SMI 1,1 ) , THH , MX , AMLER , 1 , THT AS ) 

CALL S°LCER(SM(l ,1 ),. THH, MX, AMLER, 1,DTDMS) 

CALL LININT (ZSPL IT, RSPL IT, XZ,XR,TN, 21,21, TSPLIT) 

TANS =■ RSPL IT*DTCMS 
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RIS = TSFLIT/2.0 

STGRS = 7FH1<mx)-RIS*TANS/RSPLIT-THTAS +7HHC 
AMLER=AMLER~SM(MXRL» 1 ) 

RET AS = A TAN (TANS) 

BETAS = BET AS *57 .29577 

WRITEI6.1255) AML ER , S TGR S , R I S, RO, BETAS 
4*53 WR I TE C 6 , 1 201 ) 

W’RITE(6»1239) 

WRITE (6, 12*51) 

DC 454 1 = 1, MX 

454 AE ( I ) = SMII.K ) -SM ( MX BL t K ) 

WRITEI6.1230) IAE(I),I=1,MX) 

WRITEI6.1249) 

WRITE(6,1230) ( R ( I , K >,I=l,MX) 

WRITE(6,1240) 

WRITE (6 , 1230 ) ( A C < I ) , I = l , MX ) 

WRITE(6, 1250) 

WRITF(6,1251) 

WRITE (6, 12 30) IAE( I ), 1 = 1, MX) 

WRITEI6, 1252) 

WRITE (6, 1230) (THMKl ) ,1*1, MX) 

WRITF(6»1253) 

WRITE(6*1230) <THM2( I)*.I = 1»MX) 

R I = T I ( 2 ) / 2# 

RC = TO ( 2 ) / 2 . 

STGR = TFMl(MX) 

WRITER, 1254) STGR, R I ,R0 

IF ( 2 ( MX , K ) .LT. ZSPLIT) GQ TO 456 

CALL SPLINT (Z(1,K ),SMtl,K ) , MX , Z SPL I T , 1 , AML ER ) 
CALL SPLINT (SM<1,K ),RI1,K ) , MX , AMLER , 1 , RSPL I T ) 

CALL SPLINT ISMll.K ) , THM , MX , AMLER , 1 ♦ THT A S ) 

CALL SPLC ER ( SM( 1 »K ) ,.T HM , MX , AML ER , 1 , DTDMS ) 

CALL LININT (ZSPL IT ,RSPL IT »XZ,XR» TN, 21,21,TSPLIT) 
TANS = RSPLIT*OTCMS 
RIS = TSFLIT/2.0 

STGRS = THM1(MX)-RIS*TANS/P.SPLIT-THTAS +THMC 
AMLER=AMLER-SM(MXBL,K) 

BETAS = A TAN ( TANS ) 

BETAS = EETAS*57. 29577 

WRITEI6.1255) AML ER , STGR S , R I S , RO , B ET A S 

456 WRITE(6,1202) 

WRITE(6,1239) 

WRITEI6.1251) 

DC 457 1=1, MX 

457 AE ( I ) = SMI t ,KMX )-SM(MXBL,KMX) 

WRITE (6, 1230) <AE( I ) , 1 = 1, MX) 

WRITEI6, 1249 ) 

WR ITE (6, 1230 ) ( R ! I , K MX ) , 1 = 1 , MX ) 

WR l TE ( 6 , 1 240 ) 

WRITE(6,1230) ( AC( I ) , 1 = 1, MX ) 

WRITEI6,1250) 

WRITEI6,1251) 

WRI TE(6 , 1230) (AE (I) , 1 = 1, MX ) 

WRITEI6.1252) 

WRITE I 6, 1230) I THS 1 ( I ) ,1 = l , MX ) 

WRITE(6,1253) 

WRITE (6, 1230) (THS2I I ),I = 1,MX) 

PI = T I I 2 ) / 2 . 

RC = TO ( 3 ) / 2 . 

STGR = TFS1IMX) 

WRITE (6, 12 54 ) ST GR , R I ».R0 

IF ( Z (MX ,KMX ) .LT . ZSPLIT) GO TO 459 
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call SPLINT (Z(l ,KMX),SM(1,KMX>,MX,ZSPLIT,1,AMLER) 
CALL SPLINT I SM ( 1 ,KMX ) , R ( 1 , KMX 1 , MX » AMLER , 1 , R SPL I T ) 
CALL SPLINT ( Sv ( 1 , KMX ),, THS , M X , AMLER, 1 , THTA S ) 

CALL SPLCERISMU ,KMX) ».T HS , MX » AHLER , 1 , OTOMS I 
CALL LININT (ZSPL I T , RS PL IT , XZ ,XR , TN , 2 1 , 2 1 , TSPL IT) 
TANS = R S PL I T *DT CMS 
RIS = TSFLIT/2.0 

STCRS = THS1 ( MX) -RIS*TANS/RSPLIT-THTAS +THSC 
ANLER=ANLER-SM(MXBL»KMX ) 

RETAS = ATAN(TANS) 

BETAS = BETAS *57.29577 

WRITE(6,1255) AM!L ER , STGRS , R I S, RO , B ETA S 
459 DG 460 J=l,3 
I = MXRL 


K=1 

IFU.E0.2) K = ( KMX + 1 ) / 2 
IFIJ.EO.3) K=KMX 

T1P= 1 , - ( W A ( I „K) **2+2.*W*ALM-(W*R( I ,K ) ) **2 ) /2 . /C P/TEMP 
OENSTY = T 1 P**EX PON*RHO 
C = 6.283186*R(I ,K)-XN*TT(I»K) 

WIDTH = AB(MXBL) 

I F ( J . EQ.2 ) WIDTH = AC I MX BL ) 

IF(J.E0.3) WIDTH =i AC(MXBL) 

WM = BA(2)/DENSTY/C/WIDTH 
WTFETA = ALM/RII »K)-W*RI I»K) 

RETAI(J) = ATANIWTHETA/WM) 

A A ( J ) = EETAI (J) *57.29577 
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IOIC 

I02C 

I03C 

I04C 

I05C 

ro6C 

I07C 

I08C 

I09C 

HOC 

i r ic 

112C 


CONTINUE 

WRITE (6,1170) AA 
GC TC 10 

FORMAT (4I5,6F10.4) 

FORMAT (8H0RUN NO . I 3 , 10X , 25H INPUT DATA CARD LISTING 
FORMAT (7F10.4) 

FORMAT ( 1 0X24HBCD CARDS FOR DN,WA»Z,R ) 

FORMAT ( 3 6HK STAG. SPEED OF SOUND AT INLET = ,F9.2) 
(///5X13H ITERATION NO. 13) 


FORMAT 

FORMAT 

1X5FSB 

FORMAT 

FORMAT 


(1H 6X5HALPHA9X5HRC 
SX5HSC 9X5HSD ) 
(2X10HSTREAMLINEI3) 
(9F14.6) 

Z15X5H R 


9X5HSM 


9X5HBET A 9X5HTT 


9X5HSA 


FORMAT (1HL9X5H 
114X6HTTREL ) 

FORMAT (6F19.6,F 18.6 ) 

FORMAT ( 1 8H ITERATION NO. 


15X5HWA 


15X5HPRESS14X3HWTR14X3HWL 


13, 10X , 24HMAX . STREAMLINE CHANGE = 


IF 10 . 6 ) 

1 1 3 C FORMAT (1HJ) 

116C FORMAT (12F11.4) 

117C FORMAT ( ///1HL,10X,20HINLET ANGLES - HUB ,F7.2,8H, MEANF7.2,10H 
1, SHRCLDF7.2) 

1C01 FORMAT ( 8CH 

I ) 

1007 FORMAT ( IF 0 , 3 X ,2HMX , 2X ,.3HKMX , 3X , 2HMR , 3X , 2HMZ , 6X , 1HW, 14X,2HWT, 
113X,2FXN, 12X, 3HGAM, 12X, 2HAR ) 

1008 FORMAT ( 1F0 , IX ,.4HTYPF , IX, 4H MT , 2X , 3HSR W , 1 X , 4HMXBL , 5 X , 4HTEMP , 
I11X.3HALF , 12X fr 3HRHO,12X,5HPLOSS»9X,4HANGR) 

1 C09 FORMAT ( IF 0 » 1X.4HKSTH, 1X,4HNPRT , 1X,4HITER, IX, 4H KD , 4X , 5H SFACT 
19X .6FZSPL IT, 10X, 4HBET0, 1 IX , 6HC0RFAC , 9X , 3HSSN ) 

1011 FORMAT (4I5.6G15.5) 

1027 FORMAT ( 1F0 ,4HRB =,F8.5) 

1C28 FORMAT (7G15.5) 

1029 FORMAT ( 1P0,5X*8HZS ARRAY) 

1031 FORMAT ( IPO ,5X»8HZH ARRAY) 
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1032 FORMAT ( 1F0 ,5X ,8HRS ARRAY) 

1033 FORMAT! 1F0,5X,8HRH ARRAY) 

1034 FG R M A T ( lFOf5Xf.lOHTHTA ARRAY) 

1035 FORMAT ( 1F0 ,15X,21HBLADE THICKNESS TABLE) 

1036 FORMAT! 1F0,5X,8HTN< ARRAY) 

1037 FCRMAT(1F0,5X,7HZ ARRAY) 

1038 FORMAT! 1F0,5X,.7HR ARRAY) 

1039 FORMAT! 1F0,5X,8HZT ARRAY) 

1049 FCRWATUH) 

1200 PCRVATI 1H ,20X, 3HHUB) 

1201 FORMAT! 1F1»20X»4HMEAN) 

1202 FORMAT ( 1F1*20X, 6 H SHROUD ) 

1230 FORMAT ( IF .8G16.7) 

1239 FORMAT! 1F0 , 26 HST RE AM-CH ANNEL COORDINATES) 

1240 FORMAT! 1F0,31HSTREAM-CHANNEL NORMAL THICKNESS) 

1249 FORMAT! 1F0.7HR ARRAY) 

1250 FORMAT! 1F0, 17HBL ADE COORDINATES) 

1251 FORMAT! 1F0.7HM ARRAY) 

1252 FORMAT! lF0 t 27HTHETA ARRAY BLADE SURFACE 1) 

1253 FORMAT! 1F0.27HTHETA ARRAY BLADE SURFACE 2) 

1254 FORMAT! 1F0.6HSTGR = , G 13 . 5 , 3X , 4HR I = , G 1 3. 5 , 3X, 4HR0 =,G13.5) 

1255 FORMAT! 1F0,9HSPL I TTER S , 4X , 6HML FR = , G 1 3 . 5 , 4X , 7HSTGR S =,G13.5,4X, 
14HRI =,G13.5,4X,4HR0 =i, G 1 3 . 5 , 4X , 7HBET AS =,G13.5) 

END 


tIBFTC ROOT DECK 

SUBROUTINE RUUT ! A , B , Y „X , S M f R , MX ) 

C 

C RCCT FINDS A ROOT FOR IFX-Y1 IN THE INTERVAL (A,B) 

C 

COMMON SPW 
INTEGER SRW 

DIMENSION S M ( 21 ) ,RI21) 

TCLERY = Y/50000. 

IF (SRW. EC. 21) WRITE(6, 1000) A , B » Y» TOLER Y 
XI = A 

CALL SPLINT ( SM , R f MX , X 1 , 1 . FX 1 ) 

IF (SRW. EG. 21) WRITE! 6, 1010) Xl.FXl 
X2 = P. 

10 DC 3C 1=1,15 

X = ( X 1 +> 2 ) / 2 • 

CALL SPLINT ( SM , R *MX , X , 1 , FX ) 

IF (SRW. EG. 21) WRITE(6,1010) X,FX 
IF ( (FX1-Y)*(FX-Y) .GT.O. ) GO TO 20 
X2 = X 
GO TO 30 
20 XI = X 

FX 1 = FX 
3C CONTINUE 

rF(APS(Y-FX).LT.TOLERY) RETURN 
WRITE (6,1020)A,B,Y,FX,X 
RETURN 

1C0C FORMAT (32F1INPUT ARGUMENTS FOR ROOT — A =G1 3. 5 »3X , 3HB =,G13.5, 

1 3 X ,3 HY *,G13.5,3X,8HT0LERY =, G13.5/16X, 1HX, 17X, 2HFX ) 

101C FORMAT (8X ,G16.5, 4G18.5) 

1020 FORMAT (21HRQ0T CUT OF TOL ER ANC E , 2X , 3H A = , G 16 . 5 , 2X , 3HB =,G16.5,2X, 
13HY =.,G16 .5,2X,4HFX = , G 1 6 . 5 , 2X , 3HX =,G16.5) 

END 
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0 - o ci <ji o n 


$ I BF TC SNCCTH DECK 


SUBROUTINE SMOOTH ( X , Y, XH , YH , SLOPE , SSN , INF ) 

DIMENSION X ( 2 1 > , Y ! 2 1 ) ,-XH! 21 ) » YH! 21 ) , XI ( 21 ) , Yl ( 2 1 ) » INF ( 2.1 ) , 
1SLCPE(21) 

NS = SSN 
N1=NS-1 
DC 1C I = 2 »N1 
0 = 2.0 

I F ( I • EQ . ( NS- 1 ) ) 0 = 8.0 
IF { I.eo. ( NS-2 110=4.0 
IFfl.EQ.INS— 3»» 0=2.6667 
IF (X ( 1 + 1) .EQ.XI 1-1 ) ) GO TO 5 
SLCPE1 = |Y(I+1)-Y(I-1))/(X(I+1)-X(I-1)| 

IF ( INF ( I ) .EQ. 1 ) GO TO 6 

XI m = ( ( Y( I -II -SL0PE1*X( I-in-(YH( I I -SLOPE! I J * XH €1)1 )/(SLOPe(I) 
1-SLCPE1 ) 

xi 1 1 ) = ( (xi ( n-x! 1 1 i/c i+x( 1 1 

Y 1 ( I ) = YH! I ) +SL OPE ( 1 1* ( X 1 ( I I-XHI 1 1 I 
GC TC 10 

S LOPE 1 = INFINITY 

Y 1 ( I I =SLOPE( I ) * ( X ( I-1I-XH! I ) )+YH< I I 
Y1 (I ) = ( (Y 1! I >-Y( I II/C )+Y( I 1 
X1(I)=( <X(I-1)-XCI ))/D)+X( I) 

GC TC 10 

SLOPE = INFINITY 

Y 1 { I ) = SLOPE 1M X! I ) — X { 1-1 ) I+Y! I — 1 ) 

Y1(I)=;( ( Y 1 ( I ) — Y ( I) 1/OHYI I I 
XI ( I) = X (I) 

10 CONTINUE 

DC 2C 1=2, N1 
X ( I ) = XI ( II 
20 Y! I) = Y! (I) 

RETURN 

ENC 


$ l BF TC INTCRL DECK 

SUP R CUT I N E INTGRL (X,Y,N,SUM) 

DIMENSION X ( 50 ) , Y(50) ,S( 50), A! 50 I, B! 50) ,C( 50) ,F( 50) ,W( 501 ,SB(50) , 
1 G ( 5 0 ) ,E M ( 50 1 , SUM (50) 

CCNNCN SRW 
INTEGER SRW 
DC 1C 1=2, N 
1C S ( U = X ( I ) -X ( I -1 ) 

NC = N - 1 

DC 20 1=2, NO 

All )=S( I)/6.0 

B< I )= IS ! I ) + S ( 1+1 ) ) /3.0 

C( I ) = S( H1J/6.0 

2C F( I ) = (Y( 1 + 1 )-Y( I ) )/S( 1 + 1) — (Y( I ) — Y ( I — 1) )/S( I ) 

A ( N ) = - . 5 
R( 1 ) =1.0 
8 ( N ) = 1 . 0 
cm =-.5 
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F(1)=C.0 
F(N) =0.0 
wm = en ) 

SB(1)=C(1)/W(1) 

G ( 1 ) =0 . 0 
DC 3 C 1=2, N 

W(I) = S(I)-ft(I ),*SP ( 1-1) 

SPU )=C I I) /W< I) 

3C G( I ) = {F ( i)-/s { i)*G(i-l )) /W( I ! 

EP(M=G(N > 

00 40 I =2 » N 
K=N+1-I 

4C EP(K)^C, (K)-SB(K)*EM(K + 1) 

SUP ( 1 > =C.C 
OC 50 K = 2,N 

5C SUPIK) = SUP(K-1 )+S(K )*( YIK ) +Y(K-1 ) ) / 2 . 0-S I K > ** 3* < EM I K ) +EM I K-l) ) /2 
14. C 

IFCSRVi.EC.17) WR IT £(6*1000) N , ( X ( I) , Y (II , SUM ( I) , EM (I) , I = 1 ,N ) 

RETURN 

IOCC FORMAT (17HK NO. OF POINTS =I3/10X5HX 15X5HY 15X5HSUM 
1 13X10F2ND DFR IV./I4E20.8) ) 

ENC 


tlBFTC CCNTIN DECK 

SURRCLTINE CONTIN ( WA ».WTFL , I NO * I , WT ) 

OIPFNSION SPEED! 3) ,WEIGHTC 3) 

135 GC TC (140,150,210,270,370), INO 
14 C SPEECH) = WA 

WEIGFT(l) = WTFL 
DELTA = V T / WT FL*W A-W A 

IF ! APS ( DELTA) .GT .100. ) DELTA = S IGN < 100 ., DELTA ) 
WA = CELTA+WA 

INC = 2 
RETURN 

1 5 C IF ( (WTFL-WEIGHT ( 1 ) )/(WA-SPEEDt 1 ))) 180,180,160 
16C SPEED (2 ) = WA 

DELTA = (WT-WTFL )/ (WTFL-WEIGHT ( 1 ) )*( W A- SPEEDC 1 ) ) 
IF ( APS ( CELTA ) .GT. 100. ) DELTA = S IGN ( 1 00 . , DEL T A ) 
WA = CELTA+WA 
1 66 SPEFC(l) = SPEED12) 

WEIGHT! 1 ) = WTFL 
R E T L R N 

1 7 C WRITE (6,1000) I , WTFL 
INC = 6 
RFTLRN 
16C INC = 3 

IF (WTFL.GE.WT) GO TO 140 
IF (SPEED I 1 )-WA) 190, 200, 200 

19C SPEED (2 ) = SPEED ( 1 ) 

SPEEC(l) = 2.0*SPEFD( 1 )-WA 

S p E FC ( 3 ) = WA 

WE 1CFT( 2 ) = WE IGHT ( 1 ) 

WE!GFT(3) = WTFL 
WA = SPEEC(l) 

RETURN 
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2CC SPEEC (2 ) = WA 

SPEEC ( 3 ) = SPEED ( 1 ) 

SPEED ( 1 ) = 2.0*W A-SPEED! 1 ) 

WEICHT(2> = WTFL 
WE IGHT { 3 ) = WEIGHT ( 1 ) 

WA = SPEEO(l) 

RETURN 

2 IC WE IGFT ( I ) = WTFL 

IF (WTFL.GE.WT) GO TO 140 
IF (WEIGFT(l)-WE IGHT( 2) ) 230,330,220 
22 C W E I G F T ( 3 ) * WE IGHT ( 2 ) 

WE IGHT! 2 ) = WEIGHT ( 1 ) 

SPEED (3 ) = SPEED I 2 ) 

SPEED (2 ) = SPEED ( 1 ) 

SPEED ( 1 ) = 2.0*SPEED!2)-SPEEC(3) 

WA = SPEED! 1 ) 

RETURN 

23C IF (SPFEC(3)-SPEEC<1>-10.0) 170,170,240 
24C I NO = 4 

245 IF (WFIGFTI3) -WEIGHT! 1) ) 260,260,250 
25C WA = (SPEFD!1)+SPEED(2) 1/2.0 
RETURN 

26C WA = (SPEED(3)+SPEED(2) )/2.0 
RETURN 

27C IF !SPEED !3)-SPEED! D-10.0) 170,170,280 
28C IF ( WTFL- WE IGHT! 2 ) ) 320,350,290 
29C IF (WA-SFEED! 2) ) 310,300,300 
3CC SPEED ( 1 ) - SPEED ! 2 ) 

SPEED f 2 ) = WA 
WEIGHT! 1 ) = WEIGHT ! 2 I 
WE IGHT ! 2 ) = WTFL 
GO TC 245 

3IC SPEED ! 3 ) = SPEED!2) 

S P E E D ( 2 ) => WA 

WE IGHT ( 3 ) = WEIGHT ( 2 ) 

WE IGHT ( 2 I = WTFL 
GG TC 245 

32C IF < WA-SFEED 12)) 340,330,330 
33C WE IGHT ( 3 ) = WTFL 
SPEEC t3 ) = WA 
GO TC 245 

34C WE IGHT ( 1 ) = WTFL 
SPEEC (1 ) = WA 
GC TC 245 
35C INC = 5 

IF ( WA-SFEED! 2) I 380,360,360 
36C SPEEC ( 1 ) = SPEED ( 2 ) 

WEIGHT! 1) = WEIGHT ( 2 ) 

SPEED ( 2 ) = !SPEEC!l)+SPEED(3))/2.0 
WA => SPEEDI2I 
RETURN 
37C INC = 4 

WE IGHT ( 2 ) = WTFL 
WA = (SPEEC! 1 )+SPEED(2) J/2.0 
RETURN 
38C INC = 5 

39C W F I G H T ( 3 1 = WEIGHT(2) 

S P E E C ( 3 ) = SPEED ( 2 ) 

SPEEC ( 2 ) = (SPEED!l).+ SPEED(3))/2. 

WA = SP E E 0 ( 2 ) 

RETURN 

10CC FORMAT !/12H FIXED LINE I2,12H, MAX WT = F10.6) 

ENC 
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■o 


*-IBFTC SPLCER DECK 

SUBROUTINE SPLDERI X * Y ,-N » Z , M AX , DYDX 1 

DI PENSION X(50)».Y(50)».S(50)»A(50)»B(5Q)*C(50I,F(50),W( 50 ) .S8I50), 
1GI5C ) ,EP(50),Z(50) ,DY0X( 50) 

DC 1C 1=2 * N' 

ic s< n=x( i ) —x t i-n 

N0=N- 1 

DC 20 I =2, NO 
A (I >=S(I> /6.0 
BII)afSm+S(I*l) 1/3.0 
cii >= si i + n/ 6.0 

2C F( I ) = (Y(! + 1 )-Y( I ) )/S( I + 1)-(Y( I )-Y(I-in/S( I ) 

AIN) =— . 5 
8(11=1.0 
8 I N 1 = 1 • 0 
C'( 1 ) =-• 5 
FUI=0.0 
FIN) =0.0 
W(l)=B(l) 

sB(ir=c(i i/wm 
o ( l ) =o . o 

DC 3C 1=2, N 

w( i )=e i n-A( i ) *s b i i-i i 
sbu )=c ( n/w< i) 

3c Gm = <F( h-m d*g( i-i );>/wi n 
EP(M=G(N) 

DC AC 1=2 ,N 
K=N+ 1— I 

AC EP(K ) = G(K )-SB(K) *EM(K+1) 

DC 90 1=1 , PAX 
K= 2 

tFIKIl-MlI) 60,70,70 
6C WRITE ( 6 * 1 000 ) Z I I) 

I0CC FCRPAT ( 1 7H OUT OF BLADE Z =F10.6) 

GC TC 85 

6 5 WRITE ( 6,1000 )Z(I) 

K = N 

GC TC 85 

7 C IF ( Z ( I )-X (K ) ) 85,85,80 
ec K=K+ 1 

IFIK-N) 70,70,65 

85 DYDX ( I )=~EP (K-l) *( X (K )-Z ( I ) ) **2 / 2 . 0/ S ( K ) +EM ( K ) * ( X I K- 1 ) -Z (I) )**?/2. 

10/ S (K)+< Y(K)-Y(K-1 I )/S(K)-(EP(K ) -EM (K-l ) )*$ IK 1/6.0 
8C CONTINUE 
ICC RETURN 
ENC 


5IBFTC SPLINE DECK 

SUBROUTINE SPLINE ( X , Y , N , SLO PE , EM 1 

DIMENSION X(50)*Y(50)*.SI50),A(50)*B(50)»C(50)*F(50),W(50)*SB(50), 
1G ( 50 ) ,E P ( 50 ) ,SLOPE(50) 

CCPPCN C 
INTEGER C 
DC l C I = 2 , N 
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ic s ( i ) =x ( i j -x ( i— i ) 

NC=N-1 

DC 20 1 = 2, NC 
A( I )=S( I )/6. 

B( I )= (S ( I )+S ( 1+1 ) )/3. 
c( n=S( mi/6. 

2C F(I )=(Y{ I + l)-Y( I ( )/S( 1 + 1 ) — { Y C I )-Y( 1-1) )/S( I ) 

A ( N ) =— . 5 
B ( 1 ) = 1 . 

B(N)=1. 

C(1 )=-.5 
F ( 1 ) =0 « 

F(M=0. 

W { 1 ) =B ( 1 ) 

sen )=c ( l )/w( i) 

G( 1 )=c. 

DC 30 1 = 2 * N 

MI )=P( I)-A( I )*SB< 1-1 ) 

sb < i n/wm 

3C G{I)=(F ( I )-A(I)*G( I-l))/W( n 
EF(M=G(M 
DC AO 1 = 2, A 
K = N + 1- 1 

AC EY(K)=G(K)-SB(K)*E W (K+1) 

SLCPEd )=-S<2 )/6 .*(2.*EMm+EM<2) )+(Y(2)-Y(l)I/S(2) 

DC50 1=2, N 

sc slope (i ) = sm/6. *( 2 .*em( n+Evu-m+iYm-Yi i-i ) )/sm 

IF IC.EC.13) WRITE (6,100) N , ( X ( I ) , Y ( I) , SLOPE (I) , EM (I) , I = 1 ,N ) 

ICC FCRYAT ( 2 X 1 5HNO. OF POINTS =I3/10X5HX 15X5HY 15X5HSLOPE1 5 X5H 

1EV / ( AF 20 . 8 ) ) 

return 

ENC 


. *IBfFTC LIMNT DECK 

SUBROUTINE LININT(X1,Y1,X,Y,TN,MX,MY,F) 

CCRNCN K 

DIRENSICA X(RX),Y(MY),TN(MX,MY) 

DC 1C J3=1,YX 

1C IF(X1.LE.X( J3))G0 TO 20 
J3 = RX 

2C DC 30 JA = 1 , YY 
3C IF ( Yl.LE.Yl JA ).)GO TO AO 
JA = R Y 

AC J 1 = J3-1 
J2=JA-1 

EPS1=(X1-X(J1))/(X(J3).-X(J1) ) 

EPS2=(Y1-Y(J2))/(Y(JA)-Y(J2) ) 

EPS3 = 1.-E PS1 
F p SA = 1 . -E PS 2 

F = TM J1 , J2)*EPS3*EPSA+TN(J3, J2)*EPS1*EPSA+TN( Jl, JA ) *FP S2*EPS3+' 
1TM J3,JA)*EPS1*EPS2 

IF(K.EO.IA) WRITF(6,1)X1,Y1,F,J1,J2,EPS1,EPS2 
1 FCRRAT (8H L IN I N T3F 1 0. 5 , 2 1 3 , 2F 10. 5 ) 

K = C 

RETURN 

END 
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$ I BF 1C SPLINT DECK 


SUP PCL'T I N E SPLINT < X , Y , N , Z , M AX , Y INT ) 

DIMENSION X(50),Y(50),.S(50),A(50),B<50) ,C(50),F<50) , W( 50 ) *SB l 50 ) , 
1G(5C),EN(50),Z(50>,YINT(50> 

CCFA'CN C 
INTEGER C 
DC 1C 1=2, N 

ic sm = xi n-xii-i) 

NC=N-1 

DC 20 1=2, NC 
A( I )=S(I»/6.0 
B ( I ) = { S ( I ) +S ( I+l ) )/3.0 
C(I)=S( I41J/6.0 

2C f(d = iy( h-i i-Yi i n/si i+i)-iy( n-Yt i-m/sm 

A ( N ) = -. 5 
B ! 1 ) = l • 0 
B C N ) = 1 . 0 
Cll I =- . 5 
F ( 1 ) =0 . 0 
F I N ) =C • 0 

wm = em 

SB ( 1 ) =C ( 1 ) /W 1 1 1 
Gin=o.o 

DC 30 1=2, N 

W< I ) = P< I )-A( I )*SR< 1-1 ) 

SB<!)=C(I)/Wm 

3C G(I) = (F ( IJ-A(I)*G( 1-1);) /W( I) 

DC AC 1=2 ,N 
K =N + 1 - 1 

AC EMK)=G(K )-SR(K) *EM(K + 1) 

DC PC 1 = 1 ,NAX 
K = 2 

IF(Z( I)-X I m 60,50,70 
5 C YINT(I)=V(1) 

GC TC 90 

6C IF(Z(I).LT.Il.l*Xm-.l*X(2)))WR!TE (6,1000)Z(I) 

GC TC 85 

10CC FCRBAT t 1 7 H OUT CF RANGE Z =F10.6) 

65 IF { Z ( I ) .GT. ( 1 .1*X(N)-.1*X<N-1) )) WRITE (6,1 000 ) Z ( 1 1 

K = N 

GC TC 85 

7 C IF ( Z ( I )-X (K ) ) 85,75,80 
75 YINT (I) =Y (K ) 

GC TC 90 

8 C K = K + 1 

IF(K-N) 70,70 ,.65 

P 5 YINT(I) = EP( K-l )*(X(K)-Z( I ) )**376. /SIK )+EMIK)*( ZII ) -X ( K-l I) **3/6. 
l/S(K)+ ( Y(K)/S(K)-EM(K)*S(K)/6.)*m I)-X(K-l) ) +( Y I K- 1 ) / S ( K )-EM I K-l ) 
2*S <K )/6.l*(X(K)-7(I)) 

9 C CONTINUE 

EP(N)=G(M 
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N*X A = MA>0 ( N , MAX ) 

IF(C.E0.16J WRITE (6, 1010) N, MA X » ( X ( I ) , Y (I ) , Z ( I ) , Y IN T ( I ) , I - 1 , MX A ) 
101C FCR PAT { 2 X2 IHNO. OF POINTS GIVEN =,I3,30H, NO. OF INTERPOLATED ROI 
1 NT S = « I 3 » / 10 X5HX 15X5HY 1 2X 1 1HX- IN TERPOL . 9X 1 IHY- I NTg RPOL . / ( 4 

2E2C.R) ) 

ICC RETURN 
END 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, November 16, 1971, 
132-15. 
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APPENDIX - SYMBOLS 


A coefficient, eq. (2) 

B coefficient, eq. (2) 

C coefficient, eq. (2) 

c p specific heat, J/(kg)(K); (ft)(lbf)/(slug)(°R) 

D coefficient, eq. (2) 

h enthalpy, J/kg; (ft)(lbf)/slug 

m meridional streamline distance, m; ft 

N number of blades 

2 2 

Ap” loss in relative total pressure, N/m ; lb/ft 
R gas constant, J/(kg)(K); (ft)(lbf)/(slug)(°R) 

r radius from axis of rotation, m; ft 
r c radius of curvature of a meridional streamline, m; ft 
s distance along a quasi -orthogonal, m; ft 
T temperature, K, °R 

t n blade thickness normal to mean blade shape, m; ft 
t 0 blade thickness in tangential direction, m; ft 
V absolute velocity, m/sec; ft/sec 

W relative velocity, m/sec; ft/sec 

w mass flow, kg/sec; slugs/sec 

z axial distance, m; ft 

a angle between meridional streamline and z-axis, rad 

(3 angle between relative velocity and meridional plane, rad 

y ratio of specific heats 

9 angular coordinate, rad 

2 2 

X inlet prerotation, m /sec; ft /sec 

3 3 

p density, kg/m ; slugs/ft 

xf/ angle between quasi-orthogonal and radial direction, rad 
a) rotational speed, rad/sec 


Subscripts: 

b point at which mean stream surface deviates from mean blade shape 

f flow 

h hub 

i inlet 

isen isentropic 

l leading surface 

m direction of meridional streamline 

n normal direction 

o outlet 

s shroud 

t trailing surface 

6 tangential direction 

Superscripts: 

’ absolute total conditions 

" relative total conditions 
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TABLE t. - INPUT FORM FOR QUAC 


If) 20 21 


MX KMX MR MZ 


TYPE MT SRW MXBL TEMP 


KSTH NPRT ITER KD SFACT 


THTA ARRAY 


ZT ARRAY (IF KD - 2, ZT ARRAY USED. THTA = f(ZT). MT = No. of.) 


TABLE II. - SAMPLE OUTPUT 


i 



IMPELLER 4 K E S-I&3 B 


RUN N C . 1 INPUT CATA CARD LISTING 


'KMX MP 

MZ W 

WT 

XN 

GAM 

AR 


13 21 5 

6 7P.53.98 

0.62 160E- 

-01 15.C000 

1.40000 

1716.20 


TYPE M T SRV 

MXBL TFMP 

ALM 

RHO 

PLOSS 

ANGR 


-G -0 -0 

4 51 R , 700 

-0 

0.23770E- 

-OZ 1622.80 

45.0000 


KSTF NPRT ITEP 

KD S F ACT 

ZSPLIT 

BETO 

CORFAC 

SSN 


4 5 l 

-0 2.00000 

0.12500 

-53.9500 

0.10000 

8.00000 


ZS ARRAY 

-0.754306-01 

-0.50000E-01 

-0.25000F-01 

0 

0. 99208F-02 

0.39682E-01 

0*7291 7E-01 

0.97583E-01 

0.13142 

0.16950 

0.18258 

0.18892 

0.19167 


ZH ARRAY 

-0.75430E-91 

-0.50000E-01 

-0 . 25000E-01 

0 

0.10167E-01 

0. 54000E-0 1 

0. 120P3 

C .15575 

RS ARRAY * 

0.18433 

0.201C0 

0.20475 

0.20589 

0.20617 


C. 14125 

0.14125 

0.14125 

0.14125 

0.14125 

0.14135 

0.14320 

0 • 14692 

0.15712 

0.18291 

0.20487 

0.22757 

0.25004 


RH ARRAY 

C.42375E-01 

0.42375E-01 

0.42375E-01 

0 • 42375E- 0 1 

0.42375E-01 

0.43229E-01 

0.604426-01 

0.856S5E-01 

0.12620 

0.17845 

0.20711 

0.22771 

0.25004 


THT A A R R A > 

1.C8290 

0.69770 

0.33800 

0 

-0. 12500 

-0. 50000 

-0.90000 

-I . 100CC 

-1.30000 
8L4CE THICKNESS 

-1.50000 

TABLE 

-1.60000 

-1.67500 

-1.75000 


TN ARRAY 

0.416706-02 

0.50000E-02 

0.58330E-02 

0.75000E-02 

0.75000E-02 

0.750006-02 


C.3125CE-02 

0.37500E-0? 

0.43750E-02 

0.58330E-02 

0. 75000E-02 

0.750006-02 


0.20830E-02 

0.20830E-02 

0.225006-02 

0.33330E-02 

0.500006-02 

0.583306-02 


-0 

-0 

-0 

-0 

0. 208306-02 

0.541706-02 


-c 

-0 

-0 

-0 

-0 

0.500006-02 


Z ARRAY 

-0.833C0E-03 

0 .47 250E-01 

0.94450E-01 

0.14170 

0. 18333 

0.20833 


R ARRAY 

0.4L67CE-01 

0.91 580E-01 

0.14125 

0.20833 

0.25417 





STAG. SPEEC CF SOUNC AT INLET = 1116.36 

RR = C. 2 174 7 


HERAT ICN 

NO. 

1 

MAX. 

streaml ine 

CHANGE 

= 

C. 018577 

I TERAT ICN 

*0. 

2 

MAX. 

streamline 

CHANGE 

= 

C. 015907 

ITERATION 

ISC. 

3 

MAX. 

streaml ine 

CHANGE 

= 

C. 014245 

ITERATION 

M7. 

4 

MAX. 

STREAMLINE 

CHANGE 

= 

0.012898 

ITERAT ICN 

NO. 

5 

MAX. 

streamline 

change 


C. 011650 

I TERAT ICN 

NO. 

6 

MAX. 

STRPAML INE 

CHANGE 

= 

C. 010522 

I TERATICN 

NC, 

7 

MAX. 

STREAML INE 

change 

= 

C. 009499 

ITERATICN 

NO. 

8 

MAX . 

STREAML INE 

change 

= 

C. 008590 

ITERATION 

NO. 

9 

MAX. 

STREAMLINE 

change 

= 

C. 007747 

ITERATION 

NO. 

10 

MAX. 

streamline 

change 

= 

C.0069R3 

ITERATICN 

NO. 

11 

MAX. 

streaml ine 

change 

- 

C. 006293 

ITERATICN 

NO. 

12 

MAX. 

STREAMLINE 

change 

* 

Q.00566R 

ITERATICN 

NO. 

13 

MAX. 

streaml INE 

change 

= 

C. 005104 

I TERATICN 

NO. 

14 

MAX. 

streaml INE 

change 

= 

0.004611 

ITERATICN 

NC. 

15 

MAX. 

STREAML INE 

change 

= 

C.00414B 

ITERATICN 

NO. 

16 

MAX. 

STREAMLINE 

change 

= 

C. 003731 

I TERATICN 

NC. 

17 

MAX. 

STREAML ine 

change 

= 

C. 003354 

ITERATION 

NO. 

18 

MAX. 

STREAML ine 

change 

= 

0.003015 

iteration 

NC. 

19 

MAX. 

STREAML ine 

change 

= 

C. 002709 

iteration 

NC. 

20 

MAX. 

STREAML ine 

change 

= 

0.002434 
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TABLE n. - Continued. SAMPLE OUTPUT 


ITERATION NO. 67 


ALP HA 

RC 

SM 

BETA 

TT 

SA 

SB 

SC 

SD 

STRE AML INE 1 

-0,001619 

0.001666 

0 . 

-33.192096 

0 . 

-7.071380 

8599.316772 

-0.000200 

-1731.773972 

C. 002022 

0.003332 

0.025430 

-32.066519 

0 . 

-6.649137 

8339.272705 

0.000235 

-5723.242554 

-0,006425 

-0.015122 

0.050430 

-30.721272 

0 . 

-6. 169999 

8024.158325 

-0.000692 

5073.781799 

0.C23697 

0.057244 

0.075430 

-28.558731 

0.004744 

-5.349131 

7509.221313 

0.002212 

-1434.670074 

-0.063958 

-0.357136 

0.085597 

-26.058385 

0.004834 

-4.842164 

6899.353760 

-0.005405 

-7692.848938 

4.145020 

4.201212 

0.129438 

-16.202538 

0.005290 

2.081720 

4206.509583 

-0.150864 

-1777.882187 

26.741523 

7.103331 

0.198449 

-16.851673 

0.006449 

4.734078 

3789.024353 

-2.385294 

-1276.558319 

44.955225 

7.187098 

0.241520 

-20.260381 

0.006957 

3.757438 

4238.034607 

-3.751571 

-1323.671707 

64.044206 

5.886395 

0.291126 

-25.751180 

0.006422 

1.758786 

4482.935120 

-3.613119 

-1574.753906 

79.549380 

3.979333 

0.345971 

-31.632465 

0.005817 

0.483932 

4048.283081 

-2.623682 

-518.020050 

e5. 266557 

2.981542 

0.374875 

-36.617331 

0. 006209 

0. 172996 

-2045.489014 

-2.089269 

743.802719 

88.218536 

2.065947 

0.395507 

-41.230254 

0.006235 

0.036778 

686.672424 

-1.182478 

-18.773860 

90.105220 

0.938462 

0.417838 

-53.949998 

0.005918 

-0.000597 

12901.881836 

-0.325015 

-0.086923 

STREAMLINE 6 

3.004852 

-0.147154 

0 . 

-50. 771102 

0. 

-7.610848 

12159.041016 

0.399514 

-1552.281219 

2.685177 

-0.289882 

0.025462 

-50.012261 

0 . 

-7.391540 

11959.742432 

0.346660 

-4975.400024 

1.986931 

-0.676773 

0.050483 

-48.883551 

0 . 

-7.233922 

11922.131714 

0.250962 

4686.165161 

2.757111 

1.802283 

0.075500 

-46.720125 

0.004849 

-5.564736 

11485.923950 

0.267986 

-996.231580 

4.062589 

2.746191 

0.085582 

-44.589050 

0.004847 

-4.528446 

10912.951172 

0.321631 

-6988.275818 

8.484717 

1.156943 

0.123293 

-37.185014 

0.004838 

-3.356020 

8995.205444 

0.500646 

-3405.392944 

17.853872 

5.723180 

0.175706 

-32.627673 

0.005029 

1.376789 

7146.185120 

-0.443467 

-2030.765625 

27.212084 

3.301493 

0.210372 

-29.878273 

0.005245 

0.718950 

5459.497681 

-0.369682 

-2232.488037 

44.228121 

10.548003 

0.253909 

-29.933708 

0.005329 

4.991963 

3801.527527 

-4.859242 

-2174.733856 

71.731218 

7.527054 

0.304334 

-33.680576 

0.005068 

1.557246 

3016.578857 

-4.717300 

-593.339760 

81.412612 

4.883134 

0.331991 

-36.389583 

0.005310 

0.498906 

-2376.510742 

-3.303784 

686.496490 

85.929255 

2.713154 

0.353242 

-40.335552 

0.005248 

0.112912 

574.096863 

-1.586566 

-120.49 60 5 7 

e8. 387595 

1.207301 

0.375634 

-53.949998 

0.005031 

0.011765 

12817.429688 

-0.417955 

-12.193480 


l 

R 

WA 

PRESS 

WTR 

WL 

TTREL 

streamline l 

-0.075430 

0.042375 

617.987427 

1810.882050 

628.454063 

607.520790 

527.920029 

-0.05CQ00 

0.042375 

528.864166 

1921.923386 

555.311363 

532.416969 

527.920029 

-0.025000 

0.042375 

453.919613 

2004.897369 

430.213726 

477.625500 

527.920029 

0. 

0.042375 

638.129807 

1784.119232 

646.908478 

629.351135 

527.920029 

0 .0 1C 167 

0.042375 

575.012436 

1826.039398 

612.065453 

537.959419 

527.920029 

C. 054000 

0.043229 

440.273067 

1794.442337 

449.905170 

430.640965 

528.295410 

C.12C830 

0.060442 

335.542988 

1704.843262 

389.437489 

281.640614 

537.458199 

0.155750 

0.0B5655 

419.089806 

1746.701675 

477.649139 

360.520981 

556.372025 

0.184330 

0.126200 

521.131454 

2194.479919 

643.966003 

398.296906 

600.477211 

0.201000 

0.178450 

764.574997 

3271.660889 

892.956024 

636.193970 

682.210899 

0.204750 

0.207110 

810.338303 

4495.681458 

1009.844093 

610.832512 

738.949913 

0.205890 

0.227710 

762.966583 

6003. 736694 

916.265022 

609*668144 

784.942772 

0.206170 

0.250040 

923.185478 

7195.518982 

928.774849 

917.596107 

839.720375 

STREAMLINE 6 

-0.075430 

0.079456 

811.558098 

1812.940201 

826.095238 

797.005196 

551.115196 

-0.050000 

0.080727 

748.752609 

1933.357193 

794.702576 

702.784332 

552.160957 

-0.025000 

0.081767 

698.583702 

2027.476196 

653.749702 

743.395409 

553.028755 

0. 

0.082675 

878.032661 

1729.091812 

899.370468 

856.674393 

553.794922 

0.01C065 

0.083266 

830.890221 

1782.698273 

913.205231 

748.554077 

554.299088 

0.047529 

0.087576 

704.661819 

1890.823746 

755.231674 

654.069443 

558.079300 

0.098786 

0.098518 

631.463585 

1924.358246 

700.472382 

562.426483 

568.534744 

C.13C711 

0.112028 

664.214363 

1958. 107941 

713.508934 

614.890160 

583.140831 

0.166728 

0.136486 

666.233994 

2303.128235 

763.726509 

568.733528 

614.351204 

0.192875 

0.179603 

797.847511 

3276.374908 

935.761078 

659.932220 

684.331329 

0.199133 

0.206542 

819.723518 

4457.356262 

1039.918396 

599.528473 

737.74456B 

C. 201395 

0.227673 

767.949860 

5994.851135 

924.623489 

611.276291 

784.856064 

0.202443 

0.250040 

924.466454 

7189.067261 

909 *392128 

939.4407Q4 




TABLE II. - Concluded. SAMPLE OUTPUT 


SHRCUD 

STREAN-CHAKNEl COORDINATES 
K ARRAY 


-C.7543CC0E-C1 
C . 1332545 

-0.5000000E-01 

0.1793122 

-0.2500000E-01 

0.2048316 

0 

0. 2284C03 

0. 9920799E-02 
0.2510380 

0.3968217E-01 

0.7295862E-01 

0.9791 355E-01 

R ARRAY 

C.14125C0 

0.15712CC 

0.1412500 

0.1829100 

0.1412500 

0.2048700 

0. 1412500 
0.2275700 

0. 1412500 
0.2500400 

0. 1413500 

0.1432000 

0.1469200 

STREAR-CHAANEL NCRMAL THICKNESS 
0.3221791E-C2 0.3084384E-02 

0 > 1 884080E-C 2 0.1333241E-02 

0.2948232E-02 

0.9827734E-03 

0.3039101E-02 
0. 7701051E-03 

0. 29998516-02 
0.6961409E-03 

0.3170479E-02 

0. 2734923E-02 

0.2184033E-02 

RLACE CCCRCINATES 








H ARRAY 

-C.7543000E-C1 

0.1332545 

-0. 500000 OE-Ol 
0.1793122 

-0.2 500000 E-0 1 
0.2048316 

0 

0.2284003 

0.99207996-02 

0.2510380 

0. 3968217E-01 

0.7296862E-OI 

0.9791 355E-01 

THETA ARRAY PLACE 
0 

"1 .2762058 

SURFACE 1 
0 

-1.4772978 

0 

-1.5797404 

0 

-1.6559C60 

-0.9669029E-01 

-1.7335161 

-0.4712923 

-0.8729693 

-1.0756145 

THETA ARRAY BLADE 
C 

-1. 2970669 

SURFACE 2 
0 

-1.4959751 

0 

-1.5935324 

0 

-1,6673668 

-0.1265826 

-1.7335161 

-0.5019806 

-0.9003035 

-1.0976584 

STGP = -1.73352 

R I = 0. 10415E-02 RO = 0, 

1 90280E-03 





SPLITTERS KLER 

« 0.12674 

STGRS = -0.47327 RI = 

0.11R15E-02 

RO = 0.90280E-03 

BETAS = -38 

.2841 


I MET ANGLES - HUB -30.70, MEAN -53.79, SHROUD -60.95 



TABLE III. - INPUT FOR QUAC FOR BACKSWEPT IMPELLER 


•* 

10 

15 

20 

21 30 

31 « 

41 50 

| 51 GO 

I- 

i m.K 









JACKS WEI* 

IM PE LI. HI 







MX 

KMX 

MR 

MZ 

W 

WT 

XN 

LIAM 

Alt 

13 

21 

5 

G 

7853. 98 

.06216 

15. 0 

1.40 

ITUi 2 

TYPE 

MT 

srw 

MXBL 

TEMP 

ALM 

r no 

P LOSS 

AMIR 

0 

0 

0 

4 

518. 7 

0 

. 002377 

1622. 8 

45 0 

KS I'll 

NPRT 

ITER 

KD 

SFACT 

ZSPLIT 

BE TO 

COR FAC 

SSN 

4 

5 

1 

0 

1.0 

. 208 

-53. 95 

. I 

8. 0 

ZS ARRAY 









-.07543 

-.05 


-.025 

0 

.0099208 

. 039682 

. 072917 

. 097583 

. 13142 

. 16958 

. 18258 

. 18892 

. 19167 


Zll ARRAY 









-.07543 

-.05 


-.025 

0 

.010167 

. 05400 

. 1 2083 

. 15575 

. 18433 

. 20100 

. 20475 

. 20589 

. 20617 


RS ARRAY 









. 14125 

. 14125 

. 14125 

. 14125 

. 14125 

.14135 

14 320 

. 14092 

. 15712 

. 18291 

. 20487 

. 22757 

25004 



Ri I ARRAY 


.04237.1 

. 042375 

.042375 

.042375 

.042375 

. 043229 

. 08565:i 

. 12620 

. 17845 

.20711 

. 22771 

. 25004 

THTA ARRAY 






1.0829 

.6977 

. 338 

0 

-. 125 

500 

o 

o 

-1. 300 

-1.500 

-1.600 

-1.675 

-1.750 

ZT ARRAY (IT KD 2. 

ZT ARRAY USED. THTA = f(ZT). MT -- No. of.) 









TN ARRAY 






. 004 1 6 i 

. 005000 

.005833 

. 007500 

.007500 

. 007500 

.003125 

. 003750 

.004375 

.005833 

.007500 

007500 

.002083 

. 002083 

.002250 

.003333 

. 005000 

. 005833 

0 

0 

0 

0 

. 002083 

005417 

0 

0 

0 

0 

0 

. 005000 

XZ ARRAY 






- . 000833 

.04725 

.09445 

. 1417 

. 18333 

. 20833 j 

XR ARRAY 






. 04167 

.09158 

. 14125 

20833 

. 25417 

1 
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TABLE IV. - STREAM- CHANNEL AND BLADE COORDINATES FOR BACKSWEPT IMPELLER 


NT AN 

STRFAR-CFAAAtL CCCRD1 N AT ES 


M 

ARRAY 

-C.7546577E-C1 

-0.50017C5F.-01 

-0.2500616F-C1 

0 

0. 1001300E-01 

0.4444378E-01 

0.8904854E-01 

0.119B214 


C . 16CC369 

0.2C832P3 

0.2351293 

0.2570320 

0.2794789 




R 

ARRAY 









C.10417e7 

0.1051547 

0.1058929 

0. 1064478 

0.1068066 

0.1094699 

0.1164862 

0.1256780 


C. 1442684 

0.1807066 

0.2C6C157 

0.2276400 

0.2500400 





STREAN-CFAANFL NCRMAL THICKNESS 

C.4358762E-C2 0.4259891E-02 0.4151687E-C2 0. 4 18 1409E-02 0.4095016F-02 0. 3960968E-02 0.3570975E-02 0.3193651E-02 

0.258777CE-C2 0. 1 4926 10t-02 0. 1053555E-02 0. 8437943E-03 0 . 7244833E-03 

BLACE CCCRCIMAIES 


M ARRAY 

-C *7546 5778-8 1 -0. 500 1 7C5E-0 1 -0 . 2 5006 16E-C 1 0 0. 100 1300E-0 1 0. 444437BF-01 0 . 8904B54E-0 1 0.1198214 

0 • 16CC369 0.2083283 0.2351293 0.257C320 0.2794789 


THETA ARRAY 3LACE 

SURFACE 1 







C 

0 

0 

0 

-0.8467653E-01 

-0.4606674 

-0.8628947 

-1.0642809 

-1 .2663478 

-1.4700192 

-1.5709620 

-1.6472431 

-1.7265479 




THETA ARRAY PLACE 

SURFACE 2 







0 

0 

0 

0 

-0.1290727 

-0.5030818 

-0.9008545 

-1.0994683 

-1.2974013 

-1.49373C0 

-1.5927871 

- 1 . 6665C61 

-1.7265479 




STGR = -1.72655 

R I = 0. 

14103E-02 R0 = 

0.1 5993E-02 






cn 

w 


TABLE V. - INPUT FOR QU AC FOR FLAT- VANE DIFFUSER 



10 

15 

20 

| 21 30 

[ 31 40 

| 41 50 

1 51 60 

1 61 70 

TITLF 





.. 



- - 

FLAT VANE DIFFUSER 






MX 

KMX 

MR 

MZ 

w 

WT 

XN 

GAM 

AR 

12 

21 

7 

7 

0 

.0235 

18.0 

1.667 

593. 3 

TYPE 

MT 

SRW 

MXBL 

TEMP 

ALM 

RHO 

PLOSS 

ANGR 

0 

0 

0 

1 

737.6 

77. 35 

. 0089470 

171. 4 

90.0 

KSTII 

NPRT 

ITER 

KD 

SFACT 

ZSPLIT 

BETO 

COR FAC 

SSN 

0 

5 

2 

1 

10 

.03 

41.03 

. 1 

0 

ZS ARRAY 









.00G98 

.00819 

.00900 

. 00996 

. 0105 

. 01078 

.01089 

. 01080 

.01062 

.01007 

.00882 

.00607 



ZI! ARRAY 









.02527 

.02405 

.02325 

.02228 

.02175 

.02146 

02135 

.02138 

.02163 

. 02216 

.02343 

.02617 



RS ARRAY 









. 1917 

. 1958 

. 2000 

. 2083 

. 2167 

. 2250 

. 2333 

. 2417 

. 2500 

. 2583 

. 2667 

2729 




. 1917 

. 1958 

. 2000 

. 2083 

. 2167 

. 2250 

. 2417 

. 2500 

.2583 

. 2667 

. 2729 



0 

.05166 

. 09686 

. 17174 

. 23387 

. 28606 

. 37228 

. 40823 

. 44069 

. 47054 

. 49096 



THTA ARRAY 


ZT ARRAY (IF KD 2. ZT ARRAY USED. THTA - f(ZT). MT ■- No. of.) 


. 0025 

. 0025 

. 0025 

.0025 

.0025 

.0025 

. 0025 

. 0025 

. 0025 

.0025 

.0025 

.0025 

. 0025 

0025 

. 0025 

. 0025 

.0025 

.0025 

.0025 

0025 

0025 

. 0025 

. 0025 

. 0025 

.0025 

. 0025 

0025 

0025 

. 0025 

0025 

. 0025 

.0025 

.0025 

0025 

0025 

0025 

0025 

. 002.3 

.0025 

0025 

0025 

. 0025 

0025 

0025 

0025 

0025 

. 0025 

0023 

0025 


XX AKK.U 


. 240 
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TABLE VI. - STREAM- CHANNEL AND BLADE COORDINATES FOR DIFFUSER 


Mf- AN} 

STKE AP-O AANF.L CLCHD1NATF-S 
M ARHAY 


c. 

C.5R 3CCCO--C 1 

0.4 1C0CC IE-02 
0.6 66C0C0E-01 

0.33E CCCCF-C2 
0. 74CCC00E -C 1 

0. 166CC00E-D1 
0. 812CGO0E-01 

0.25UCO0OF-01 

0. 3330000F-01 

0.4160000E-01 

0.5000000F-01 

R ARRAY 

C . 191 7CGC 
C.2.5CCCCC 

0. 1 95RRCC) 
0.25P30CO 

0. 2C0CCCC 
0. 2667CC0 

0.2CR3C00 
C. 2729C00 

0.2157000 

0. 2250C00 

0.2333000 

0.2417000 

STREAR-CFAI'M.'L NC R R A1 THICKNESS 
C.911K, 576H-C3 C . 1 94O0C 3E-03 

0 .5513 1531- -C 3 0. ft 04 896 Bt-03 

0.7129374C-C3 
0. 7394151C-C3 

0.6164 340E-03 
0. 9687544E-03 

0.6628559E-03 

0.5343616E-03 

0.5232948E-03 

0.5 265 550E-03 

BL AC e CECPCI'IAIFS 








M ARRAY 
C 

C.583CCCCE-C1 

0.4 1C0DC 1F-02 
0.666CCCOE-01 

0.8300000C-C2 

0.75CC0Q0E-C1 

0. 166CC0OE-O1 
0.B12CC00E-01 

O.26OCO00E-01 

0. 3330000E-01 

0.4160000E-01 

0.5000000E-01 

THETA ARRAY hlacf 
C 

C. 39(17394 

SURFACE 1 
0. 5 C94 88 8F- 0 1 
0.4 30746 1 

0.G419942E-C1 

0.4601997 

0. 166F712 
0.4702807 

0.2274512 

0. 2785887 

0.3233202 

0.3633232 

THETA ARRAY HLACE 
0 

C.3844C31 

SURFACE 2 
0. 19C5362E-01 
0.4 17 3164 

0.66203CPF-C1 

0.4476628 

0.1432913 

0.47G2807 

0.2069713 

0.2602138 

0.3065823 

0.3479192 

STGR = C. 47026 

R I = 0. 1250CF-C2 R0 = 0. 

1 2 600E-02 







TABLE VII. - INPUT FOR QUAC FOR RADIAL IMPELLER 


5 1 I u | 1 "1 20 1 21 30|31 40 1 41 5ojT l 00 1 PI 7(7 


TITLE 



HAiMAJ IMPELLER 

MX 

KMX 

MR 

MZ 

W 

WT 

XN 

GAM 

AR 

10 

21 

G 

10 

4031.7 

.018975 

15. 0 

1. 667 

1246. 1 

TYPE 

MT 

SRW 

MXBL 

TEMP 

ALM 

RHO 

PI OSS 

ANtiR 

0 

11 

0 

1 

536. 0 

0 

. 001293 

164. 0 

45.0 

KSTH 

N PUT 

ITER 

KD 

SFACT 

ZSPLIT 

BETO 

COR FAC 

SSN 

0 

5 

1 

2 

1.0 

. 16 

-21.0 

. 1 

0 

ZS ARRAY 

1 0 


.02225 

.043083 

. 05975 

. 076416 

. 088916 

. 10558 

I .11907 

. 13108 

. 13317 





Zll ARRAY 

1 0 


. 02225 

.051416 

.076416 

. 10142 

. 11967 

. 1 3642 

| . 14042 

. 15142 

. 15142 





RS ARRAY 

. 14592 

. 14600 

. 14700 

. 14975 

. 15475 

. 16067 

. 17325 

. 19107 

.2250 

. 2490 





RH ARRAY 

.077410 

.0775 

.0805 

.088166 

. 1035 

. 12217 

. 15033 

. 17892 

. 22283 

.2490 






THTA ARRAY 

I 0 


-.078 

1785 

-.2659 

-.3745 

-.4802 

-.5519 

-.5924 

-.6069 

-.6120 

-.6120 




ZT ARRAY (IF KD 2. 

ZT ARRAY USED. THTA - f(ZT). MT -- No. of. ) 




1 0 


.0055833 

.013917 

.02225 

.03475 

.051416 

.068083 

. 08475 

.09725 

.11433 

. 15142 




TN ARRAY 

. 002667 

. 003417 

.00475 

. 009083 

.01075 

.01167 

.0125 

.0125 

.0125 

.0125 





. 002083 

. 002500 

.003458 

.00675 

. 00825 

.00925 

.01017 

.01083 

.01146 

.01146 





.00125 

.001417 

.001833 

.003792 

. 005167 

. 006208 

.007167 

.0080 

.008708 

.01012 





1 0.0 


0.0 


0.0 

.001417 

.002708 

. 003792 

. 004833 

. 00575 

. 00650 

.007833 





0.0 


0.0 


0.0 

0.0 

0.0 

. 001375 

. 00250 

.00350 

. 00425 

. 005667 





0.0 


0.0 


0.0 

0.0 

0.0 

0. 0 

0.0 

.001107 

.001958 

. 003208 





XZ ARRAY 

-.000833 1 

.005583 

.018083 

. 05975 

. 08475 

. 10142 

. 11433 

. 123 

i! 

. 13317 

. 15167 





XR ARRAY 

.075 

. 10833 

. 1500 

. 18333 

. 21667 

. 251fiR 
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